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Abstract  of  Thesis  Presented  to  the  Graduate  School 
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Requirements  for  the  Degree  of  Master  of  Arts 

ABSOLUTE  SENSITIVITY  AND  FREQUENCY  RESOLUTION 
AT  HIGH  AUDITORY  FREQUENCIES 

By 

Bin  Zhou 
August  1994 

Chairperson:  Dr.  David  M.  Green 
Major  Department:  Psychology 

The  present  studies  focus  on  the  auditory  perception  of 
sound  at  frequencies  above  8 kHz.  In  the  first  study,  an 
audiometry  system  was  designed  to  investigate  the  effect  of 
standing  waves  on  estimates  of  the  high-frequency  sensitivity 
of  the  ear.  In  the  second  study,  the  intra-subject 
reliability  of  high-frequency  thresholds  was  measured  for  two 
types  of  sound  sources.  In  the  third  study,  the  frequency 
resolution  of  the  auditory  system  was  estimated  for 
frequencies  above  10  kHz  using  two  different  methods. 

The  purpose  of  the  first  study  was  to  determine  the 
effect  of  standing  waves  in  the  ear  canal  on  absolute 
thresholds.  A standing  wave  field  was  created  by  inserting  a 
hard-wall  material  into  a long  canal-size  plastic  tube  which 
conveyed  sound  from  the  transducer  to  the  ear  canal . The 
results  showed  that,  in  the  frequency  range  from  5 to  16  kHz, 


v 


the  change  in  absolute  threshold  as  a function  of  frequency 
caused  by  the  standing  waves  was  about  ± 5.0  dB . The 
predicted  sound  pressure  at  the  eardrum  correlated  well  with 
the  changing  pattern  of  the  absolute  threshold  at  frequencies 
below  13  kHz. 

In  the  second  study,  the  variation  of  intra-subject 
thresholds  caused  by  transducer  placement  was  estimated  for 
two  types  of  transducers,  the  hard-wall  source  used  in  the 
first  study  and  several  conventional  headphones.  As  expected, 
it  was  found  that  the  variation  in  absolute  threshold  was  much 
larger  for  signals  above  7 kHz  than  for  a 1-kHz  signal.  The 
results  also  showed  that  the  variation  caused  by  transducer 
placement  was  larger  for  the  hard-wall  source  than  for  the 
conventional  headphones. 

In  the  third  study,  two  methods  were  used  to  estimate  the 
frequency  resolution  for  frequencies  above  10  kHz.  The  first 
method  employed  notched-noise  maskers.  The  results  showed 
that  the  equivalent  rectangular  bandwidths  of  the  estimated 
auditory  filters  at  12  and  14  kHz  were  close  to  the 
extrapolated  values  from  previous  studies.  The  second  method 
was  based  on  the  monaural  phase  discrimination  of 
amplitude-modulated  and  quasi -frequency-modulated  signals. 
The  obtained  critical  bandwidths  at  frequencies  above  8 kHz 
were  smaller  than  the  critical  bandwidths  obtained  with  the 
masking  method. 


vi 


Chapter  1 
INTRODUCTION 

The  present  studies  focus  on  the  auditory  perception  of 
sound  at  frequencies  above  8000  Hz.  An  audiometry  system  was 
designed  to  investigate  the  effect  of  standing  waves  on  the 
estimates  of  high-frequency  sensitivity  of  the  ear.  The 
intra-subject  reliability  of  high-frequency  threshold 
estimates  was  measured  for  two  types  of  sound  sources.  The 
frequency  resolution  of  the  auditory  system  was  also  estimated 
for  frequencies  above  10  kHz. 

Traditional  audiometry  tests  cover  only  the  frequency 
range  from  125  to  8000  Hz,  although  the  human  cochlea  responds 
to  frequencies  as  high  as  24000  Hz  (Wegel,  1932;  Rosen, 
Plester,  El-Moffy,  & Rosen,  1964)  . There  are  few  data  on 
hearing  above  8000  Hz.  One  reason  for  the  neglect  is  that 
much  auditory  communication  depends  only  on  lower  frequencies. 
For  example,  speech  can  be  easily  understood  over  a telephone 
system  that  does  not  reproduce  frequencies  above  3000  Hz.  In 
recent  years,  however,  there  has  been  an  increased  interest  in 
the  measurement  of  high-frequency  hearing. 

Studies  have  shown  that  the  hearing  loss  due  to  ototoxic 
drug  use  or  intensive  noise  exposure  was  first  apparent  and 
most  severe  in  the  high-frequency  region  (Jacobson,  Downs,  & 
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Fletcher,  1969;  Fausti,  Schechter,  Rappaport,  Frey,  & Mass, 
1984;  van  Zeijl,  Conijn,  Rodenburg,  Tange,  & Brocaar,  1984; 
Dreschler,  van  der  Hulst,  Tange,  & Urbanus,  1985) . The  high- 
frequency  threshold  measurement  could  be  used  as  a good  early 
warning  system  for  various  hearing  losses.  Others  also 
suggested  that  the  measurement  of  frequency  resolution  and 
temporal  integration  at  high  frequencies  should  also  be 
included  during  administration  of  ototoxic  drugs  (Collins, 
Cullen,  & Berlin,  1981;  Henry,  Fast,  Nguyen,  Paolinelli,  & 
Ayars,  1985;  Viemeister,  1982;  Viemeister  & Bacon,  1983). 

Despite  these  diverse  uses,  there  is  no  standardized 
system  for  the  calibration  and  measurement  of  high-frequency 
thresholds.  Three  most  popular  high-frequency  audiometers 
were  designed  by  Rudmose  (Fletcher,  1965),  Fausti,  Frey, 
Erickson,  Rappaport,  and  Cleary  (1979),  and  Stevens, 
Berkovitz,  Kidd,  and  Green  (1987)  . The  lack  of 
standardization  for  these  systems  is  manifested  by  the 
different  methods  of  delivering  sound  to  the  human  ear  canal. 
A metal  tube,  supraaural  headphones  and  a long  plastic  tube 
were  used  as  the  means  of  conveying  sound  to  the  listener. 

The  essential  problem  with  the  measurement  of  high- 
frequency  hearing  arises  from  the  fact  that  standing  waves  are 
generated  in  the  ear  canal  for  frequencies  above  3500  Hz. 
These  give  rise  to  resonances  when  the  ear  is  driven  by  a 
conventional  type  of  transducer,  making  calibration  difficult. 
Further,  because  of  the  small  wavelength  of  sound  at  high 
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frequencies  and  the  slope  of  the  eardrum  in  the  ear  canal, 
there  is  a distribution  of  sound  pressure  over  the  surface  of 
the  membrane,  rather  than  a single  sound  pressure.  There  is 
no  consensus  as  to  how  this  distribution  can  best  be  related 
to  the  audiometric  threshold. 

One  goal  of  the  present  study  was  to  investigate  how 
standing  waves  in  the  ear  canal  affect  high-frequency  pure- 
tone  thresholds.  The  experiments  presented  in  Chapter  2 were 
designed  to  address  the  issue.  Two  sources  were  used  to 
estimate  absolute  thresholds  measured  at  frequencies  between 
5 and  16  kHz.  One  source  was  a long  lossy  tube  that  produced 
essentially  no  reflections  at  the  entrance  to  the  ear  canal. 
The  second  source  was  a hard-wall  material  that  produced 
strong  reflections  at  the  entrance  to  the  ear  canal.  The 
difference  in  threshold  obtained  with  the  two  different 
driving  sources  revealed  the  effect  of  standing  waves  in  the 
ear  canal . These  standing  waves  produced  large  changes  in 
sound-pressure  level  as  a function  of  frequency  at  the 
eardrum.  The  pressure  measured  at  the  entrance  of  the  ear 
canal  was  used  to  predict  the  pressure  at  the  eardrum.  This 
calculated  eardrum  pressure  was  compared  with  the  behavioral 
thresholds . 

Because  standing  waves  result  in  different  sound-pressure 
levels  at  the  eardrum  for  different  positioning  of  the 
transducer,  the  different  sound-pressure  levels  produces 
variations  between  the  test-retest  results.  Several  factors 
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that  influence  the  intra-subject  reliability  of  high-frequency 
threshold  estimates  were  investigated  in  Chapter  3 . The 
reliability  of  high-frequency  audiometry  is  determined  by  two 
variances  (Stelmachowicz , Beauchaine,  Kalberer,  Kelly,  & 
Jesteadt,  1989),  threshold  variance  and  fitting  variance. 
Threshold  variance  is  determined  largely  by  the  number  of 
trials  tested,  the  slope  of  the  underlying  psychometric 
function,  the  stability  of  the  observer,  and  the 
psychophysical  procedure.  Fitting  variance  is  closely 
associated  with  the  existence  of  standing  waves  in  the  ear 
canal.  It  is  determined  by  the  frequency  of  the  sound  source 
and  the  amount  by  which  the  sound  source  moves  during  the 
test.  In  two  experiments,  the  intra- subject  threshold 
variance  and  fitting  variance  were  evaluated  for  two  types  of 
sound  sources,  the  hard-wall  source  and  the  headphones.  In 
another  experiment,  the  effect  of  the  slope  of  the 
psychometric  function  on  threshold  variance  was  investigated. 

In  Chapter  4,  we  estimated  the  frequency  resolution  of 
the  ear  for  frequencies  above  10  kHz.  Measurements  of 
frequency  resolution  have  been  undertaken  by  many  researchers. 
No  data,  however,  have  been  obtained  for  frequencies  above  10 
kHz.  One  of  the  problems  is  that  it  is  difficult  to  generate 
these  high-frequency  stimuli  with  controlled  spectra  using 
traditional  headphones.  Another  problem  is  the  variability  in 
the  high-frequency  threshold  estimates  introduced  by  the 
fitting  variance  using  traditional  headphones.  These  two 


5 


problems  can  be  avoided  by  using  the  audiometry  system 
designed  in  the  present  study.  Two  methods  were  employed  to 
estimate  the  frequency  resolution  of  the  ear.  The  first 
method  employs  notched-noise  maskers  as  described  by  Patterson 
(1976)  and  Patterson  and  Nimmo-Smith  (1980)  . The  second 
method  is  based  on  the  monaural  phase  discrimination  of 
amplitude-modulated  and  quasi -frequency-modulated  signals 
(Goldstein,  1967) . 


Chapter  2 

EFFECT  OF  STANDING  WAVES  ON  HIGH-FREQUENCY  THRESHOLDS 

Introduction 

The  clinical  significance  of  high-frequency  audiometry 
became  apparent  as  research  demonstrated  that  many  forms  of 
cochlear  damage,  including  noise  exposure,  diseases  and 
ototoxicity,  were  first  evident  at  the  basal  end  of  the 
cochlea  (Fletcher,  Cairns,  Collins,  &Endicott,  1967;  Jacobson 
et  al . , 1969;  Fausti,  Erickson,  Frey,  Rappaport , & Schechter, 
1981).  Fausti  et  al . (1981)  observed  that  excessive  noise 

exposure  caused  the  most  severe  hearing  loss  in  the  high- 
frequency  region.  Fletcher  et  al . (1967)  tested  26  patients 

with  meningitis.  Compared  to  a control  group,  the  meningitis 
patients  had  greater  hearing  loss  above  6 kHz,  but  no 
statistically  significant  deterioration  of  hearing  below  6 
kHz.  Ototoxic  cochlear  damage  as  observed  in  animals  showed 
most  hair  cell  loss  originated  at  the  basal  turn  of  the 
cochlea,  which  encodes  higher  frequencies,  and  progressed 
toward  the  apical  end,  which  encodes  lower  frequencies  (Tange, 
Conijn,  van  Zeijl,  & Huizing,  1982a) . Animal  audiograms 
showed  that  hearing  loss  induced  by  cis-diammine- 
dichloroplatinum- II  (DDP)  was  more  severe  at  high  frequencies 
than  at  low  frequencies  (Tange,  Conijn,  & van  Zeijl,  1982b) . 
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For  human  subjects,  the  audiogram  before  and  after  treatment 
with  DDP  showed  no  difference  in  the  frequency  range  of  250  to 
4000  Hz.  For  8000  Hz,  however,  over  40%  of  the  patients 
showed  more  than  a 20-dB  hearing  loss  (van  Zeijl  et  al . , 
1984)  . Studies  of  hearing  sensitivity  above  8000  Hz  both 
before  and  after  treatment  with  different  ototoxic  drugs  have 
shown  high-frequency  audiometry  to  be  a good  early-warning 
system  for  ototoxicity  (Jacobson  et  al . , 1969;  Dreschler  et 
al . , 1985;  Fausti  et  al . , 1984). 

Despite  the  clinical  significance  of  high-frequency 
audiometry,  there  is,  unfortunately,  no  standard  method  of 
testing  and  calibrating  high-frequency  thresholds.  Over  the 
low-frequency  range  of  standard  audiometric  testing  (i.e.,  125 
Hz  to  2000  Hz) , the  wavelength  of  sound  ranges  from  2774  mm  to 
172  mm.  These  wavelengths  are  much  greater  than  the  length  of 
the  average  ear  canal,  which  is  about  26  mm.  Because  the 
wavelength  is  so  long,  the  interaction  of  incident  and 
reflected  waves  produces  little  change  in  sound  pressure  over 
the  length  of  the  ear  canal.  At  higher  frequencies,  the 
interaction  of  incident  and  reflected  waves  creates  points  of 
maximum  and  minimum  pressure  along  the  ear  canal  called 
standing  waves  (Stinson,  1990) . Standing  waves  are  apparent 
when  the  quarter  wavelength  of  the  signal  is  smaller  than  the 
length  of  the  ear  canal,  which  for  the  average  human  ear  means 
frequencies  greater  than  about  3500  Hz.  The  existence  of 
standing  waves  also  implies  that  for  a fixed  input  level  the 
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sound  pressure  will  change  as  a function  of  frequency  at  any 
fixed  location  in  the  ear  canal. 

The  first  problem  encountered  in  measuring  high-frequency 
thresholds  is  how  to  deliver  sound  to  the  ear  canal.  Most 
traditional  audiometry  systems  use  supraaural  or  circumaural 
headphones  as  transducers.  Accurate  estimates  of  pure-tone 
thresholds  would  require  the  same  sound  pressure  in  the  ear 
canal  over  repeated  measurements.  The  pressure  in  the  ear 
canal  depends  upon  the  position  of  the  transducer.  For 
example,  Shaw  (1966)  demonstrated  that  three  successive 
transducer  fittings  (TDH-39)  resulted  in  pressure  changes  in 
the  ear  canal  of  7 dB  at  10  kHz.  Even  when  the  transducer 
placement  and  the  input  level  are  fixed,  the  presence  of 
standing  waves  means  that  eardrum  sound  pressure  will  change 
as  a function  of  frequency. 

The  second  problem  involves  how  the  sound-pressure  level 
should  be  calibrated.  In  most  cases,  this  is  equivalent  to 
the  question  of  where  to  place  the  calibration  microphone.  At 
low  frequencies,  the  measurement  location  does  not  matter, 
because  the  sound  pressure  level  is  approximately  uniform  over 
the  entire  ear  canal.  At  high  frequencies,  however,  different 
measurement  locations  result  in  different  calibration  values, 
because  of  the  change  in  pressure  associated  with  the  standing 
waves.  An  additional  calibration  problem  is  the  geometry  of 
the  tympanic  membrane.  The  tympanic  membrane  terminates  the 
ear  canal  at  an  angle  of  approximately  30  degrees  and  has  a 
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length  of  8-10  mm  along  the  canal  (Johansen,  1975)  . At  15 
kHz,  the  wavelength  is  22  mm,  so  the  change  in  pressure  over 
the  8-10  mm  eardrum  ranges  from  maximum  to  minimum  pressure. 
Using  a replica  of  the  ear  canal,  Stinson  (1985)  found  that 
the  sound  pressure  over  the  surface  of  the  eardrum  varied  by 
as  much  as  20  dB  at  a frequency  of  15  kHz.  Khanna  and  Stinson 
(1985)  also  observed  large  changes  of  sound  pressure  over  the 
surface  of  the  eardrum  in  the  cat.  These  results  imply  that 
at  high  frequencies  it  is  difficult  to  define  or  to  measure  a 
single  sound  pressure  at  the  eardrum. 

This  study  investigates  the  effect  of  standing  waves  on 
absolute  thresholds  measured  as  a function  of  frequency, 
especially  at  frequencies  above  8 kHz.  Before  introducing  the 
high-frequency  audiometry  system  developed  for  this  study, 
several  high-frequency  audiometry  systems  will  be  reviewed  to 
examine  how  these  systems  deal  with  the  problems  caused  by 
standing  waves  in  the  ear  canal . 

The  most  popular  high-frequency  audiometer  in  the  1960s 
and  1970s,  was  designed  by  Rudmose  (Fletcher,  1965)  . It 
covered  the  frequency  range  from  4000  to  18,000  Hz.  It  was 
widely  used  and  evaluated  by  other  researchers  (Fletcher, 
1965;  Zislis  & Fletcher,  1966;  Northern,  Downs,  Rudmose, 
Glorig,  & Fletcher,  1972;  Harris  & Myers,  1971)  . The  Rudmose 
audiometer  used  a Bekesy-type  microphone  as  the  transducer. 
Sound  was  coupled  to  the  ear  canal  through  a metal  tube  with 
a plastic  tip.  Because  the  transducer  was  hand  held,  the 
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sound  source  was  placed  at  various  distances  from  the  eardrum, 
and  the  resulting  test  data  were  widely  variable.  Another 
problem  was  the  calibration  procedure.  In  early  experiments 
with  this  device  (Fletcher,  1965) , sound  pressure  was 
calibrated  in  the  open  air  with  a 1/2  inch  Briiel  and  Kjaer 
microphone.  In  a later  study,  Northern  et  al . (1972) 
positioned  a Briiel  and  Kjaer  1/8  inch  probe  microphone  close 
to  the  eardrum  to  obtain  the  "eardrum  threshold  SPL."  The 
existence  of  the  probe  microphone  influenced  the  sound  field 
in  the  ear  canal,  and  different  positions  of  the  probe 
microphone  led  to  different  calibration  results. 

Later,  Fausti  et  al . (1979)  designed  another  high- 
frequency  audiometry  system.  A supraaural  headphone  (Koss 
HV/lA)  was  used  as  the  transducer  and  a flat -plate  coupler  was 
designed  for  the  calibration  of  the  sound  pressure.  As 
discussed  above,  the  use  of  a supraaural  headphone  will 
generally  lead  to  standing  wave  patterns  in  the  ear  canal 
that,  in  turn,  will  influence  the  measurement  of  high- 
frequency  thresholds . 

Another  high-frequency  audiometry  system  was  the  quasi - 
free-field  system  developed  by  Osterhammel,  Osterhammel,  and 
Terkildsen  (1977)  . It  has  been  used  by  several  researchers  to 
test  high-frequency  thresholds  (De  Seta,  Bertoli,  & Filipo, 
1985;  Henry  et  al . , 1985)  . In  this  system,  sound  delivered  to 
the  ear  canal  was  calibrated  at  the  entrance  to  the  canal 


rather  than  the  eardrum. 
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The  only  high-frequency  audiometry  system  which 
incorporated  differences  in  ear-canal  size  and  geometry  was 
designed  by  Stevens  et  al . (1987) . In  this  system,  sound  was 
delivered  to  the  ear  canal  by  a long  plastic  tube  yielding  a 
source  impedance  that  was  approximately  matched  to  the 
characteristic  impedance  of  the  ear  canal.  This  approach 
specified  threshold  measurements  in  dB  SPL  at  the  tympanic 
membrane  by  calibrating  the  sound  pressure  according  to  the 
standing  wave  patterns  in  the  ear  canal.  However,  tests 
showed  that  the  reliability  of  thresholds  obtained  with  this 
system  was  not  superior  to  the  Koss  headphone  system,  due  to 
the  variance  associated  with  earpiece  replacement  and  the 
calibration  process  (Stelmachowicz  et  al . , 1989). 

At  present,  most  practical  high-frequency  audiometry 
systems  still  use  supraaural  or  circumaural  headphones.  It  is 
interesting  to  note  that,  despite  the  presence  of  standing 
waves,  the  reliability  of  high-frequency  thresholds  is 
reported  to  be  as  good  as  that  obtained  with  low-frequency 
thresholds.  For  example,  Fletcher  (1965),  using  the  Rudmose 
audiometer,  and  Dreisbach  and  Frank  (1991) , using  a Beltone 
2000  audiometer,  found  that  intra-subj ect  high-frequency 
thresholds  were  as  reliable  as  low-frequency  thresholds. 
Testing  more  than  100  patients  treated  with  ototoxic  drugs, 
Dreschler  et  al . (1985)  also  reported  that  the  intra-subject 
standard  deviations  for  high-frequency  threshold  measurements 
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were  nearly  as  good  as  those  obtained  from  the  low-frequency 
measurements . 

Our  primary  interest  is  to  determine  how  standing  waves 
in  the  ear  canal  influence  high-frequency  thresholds.  The 
procedure  is  to  fix  the  distance  between  the  sound  source  and 
the  eardrum  and  to  change  the  frequency  of  the  input  signal. 
Even  with  the  same  input  signal  level,  the  sound  pressure  at 
the  eardrum  will  vary  as  a function  of  frequency.  This 
variation  in  the  sound-pressure  level  should  exert  a 
discernible  effect  on  absolute  thresholds. 

In  the  remainder  of  the  chapter,  a system  is  described 
that  allows  quantitative  prediction  of  the  sound  pressure  at 
the  eardrum  from  measurements  made  at  the  entrance  to  the  ear 
canal.  Next,  the  absolute  thresholds  for  the  listener  are 
measured  with  two  types  of  sound  sources.  Finally,  the 
predicted  change  in  sound  pressure  at  the  eardrum  is  compared 
with  the  estimate  of  threshold  measured  as  a function  of 
frequency. 

Estimation  of  Sound  Pressure  in  the  Ear  Canal 

Two  different  sources  were  used  to  deliver  sound  to  the 
ear  canal.  One  source  was  a long  lossy  tube.  The  other 
source  was  a hard-wall  material.  The  ear  canal  was  treated  as 
a simple  tube,  and  equations  that  related  sound  pressures  at 
various  locations  in  the  tube  were  derived.  The  equations 
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were  used  to  predict  the  pressure  at  the  eardrum  from  the 
measured  pressure  at  the  entrance  to  the  ear  canal. 

Experimental  Apparatus 

Figure  2-1  shows  a block  diagram  of  the  experimental 
apparatus.  Stimuli  were  generated  digitally  and  were 
delivered  via  a digital -to-analog  converter  (TDT  QDA2 ) to  a 
power  amplifier  (Crown  D60)  . The  power  amplifier  drove  a 
transducer  (Radio  Shack  40-1310B)  , which  was  coupled  to  the 
ear  with  a fiberglass  horn  and  a long  lossy  plastic  tube  (2  m 
long,  8 mm  diameter)  . The  tube  was  attached  to  a custom 
earmold,  which  was  placed  in  the  listener's  ear.  A small 
microphone  (Knowles  Electronics  EA-1954)  was  mounted  about  7 
cm  from  the  earmold  in  the  plastic  tube.  The  microphone 
output  was  amplified,  digitally  converted  (TDT  QAD1) , and 
stored  in  the  computer. 

Sound  Delivery  Systems 

Tube-earmold-canal  system 

A custom  earmold  was  constructed  for  each  subject  that 
allowed  precise  coupling  of  the  short  plastic  tube  to  the  ear 
canal  and  prevented  sound  leakage  at  that  juncture.  A short 
plastic  tube  (=  7 cm  long,  8 mm  diameter)  was  permanently 
attached  to  the  earmold.  The  tube  diameter  was  chosen  to  be 
equal  to  the  size  of  an  average  ear  canal.  The  earmold  had  a 
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Figure  2-1. 


Block  diagram  of  the  experimental  apparatus 
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hole  cut  in  its  middle  that  accepted  the  short  plastic  tube. 
The  earmold  was  formed  so  that,  when  in  place,  the  plastic 
tube  and  ear  canal  comprised  a tube  of  nearly  uniform  cross- 
sectional  area.  This  tube-earmold-canal  system  was  driven  by 
one  of  the  two  sound  delivery  systems. 

Long  lossy  tube 

In  one  set  of  measurements,  a long  (2  meter)  lossy  tube 
was  connected  to  the  tube-earmold-canal  system.  The  diameter 
of  this  tube  was  equal  to  the  diameter  of  the  short  tube 
connected  to  the  earmold.  The  measured  loss  in  the  tube  at  5 
kHz  was  about  12  dB  per  meter.  Thus,  the  pressure  level  of 
the  reflected  sound  from  the  transducer  was  at  least  48  dB 
lower  than  that  of  the  direct  sound  pressure  level.  At  the 
eardrum,  there  was  essentially  no  change  in  pressure  as  a 
function  of  frequency  for  a constant  input  level. 

Hard-wall  source 

In  another  set  of  measurements,  an  aluminum  cylinder  was 
inserted  at  the  end  of  the  long  lossy  tube  to  produce  the 
hard-wall  source.  The  aluminum  cylinder  was  located  about  7 
cm  from  the  custom  earmold  or  about  9.7  cm  from  the  eardrum 
(see  Fig.  2-2a)  . There  was  a very  small  hole  in  the  middle  of 
the  cylinder  that  allowed  sound  to  pass  through  with  an 
attenuation  of  about  15  dB . Because  the  hard  wall  of  the 
cylinder  reflected  sound  back  into  the  ear  canal,  substantial 
standing  waves  occurred  with  corresponding  resonances  and 
anti-resonances  as  a function  of  frequency.  Thus,  at  the 
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Figure  2-2.  (a)  schematic  of  the  area  between  the  hard-wall 
material  and  the  eardrum  --  tube-earmold-canal  system,  (b) 
diagram  of  the  physical  model  of  the  tube-earmold-canal 
system . 
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eardrum,  there  would  be  considerable  change  in  sound  pressure 
as  a function  of  frequency  even  if  the  input  sound  pressure 
was  constant . 

Physical  Model  of  the  Tube-earmold-canal  System 

A good  first-order  model  of  the  entire  tube-earmold-canal 
system  is  a uniform  tube  of  a constant  cross-sectional  area 
(see  Fig.  2-2b)  . The  length  of  the  model  tube  is  9.7  cm, 
which  approximated  the  distance  from  the  hard-wall  material  to 
the  eardrum  of  the  listener.  The  position  near  the  right  end 
of  the  tube  is  called  the  'eardrum'  position,  and  the  opposite 
end  is  called  the  'entrance'  position.  The  eardrum  position 
is  terminated  with  an  aluminum  surface  with  an  impedance 
denoted  by  Ze. 

Using  Golav  Codes  to  Measure  the  Frequency  Response 

Before  discussing  the  frequency  response  of  the  physical 
model  and  the  human  ear  canal,  we  will  first  describe  the 
signal  used  to  obtain  the  frequency  response.  Measurements  of 
frequency  responses  with,  literally,  a brief  impulse  are 
complicated  by  the  high  peak  factor  of  the  stimulus.  That  is, 
the  small  ratio  of  RMS-to-peak  pressure  can  result  in  a poor 
signal-to-noise  ratio.  This,  in  turn,  forces  the  experimenter 
to  average  responses  over  a large  number  of  stimulus 
presentations.  There  have  been  several  efforts  to  devise  a 
broadband  stimulus  that  minimize  the  peak  factor  (Schroeder, 
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1970;  Pumplin,  1985)  . We  have  adopted  a type  of  digital  probe 
signal,  the  Golay  codes,  described  by  Golay  (1961)  and 
suggested  for  use  in  impulse  response  measurements  by  Foster 
(1986) . Golay  codes  consist  of  pairs  of  signals,  each 
comprising  a series  of  positive-  and  negative-going  pulses. 
We  will  demonstrate  that  a measurement  made  with  a single 
presentation  of  a pair  of  Golay  codes  of  length  L can  approach 
the  theoretical  improvement  in  signal -to-noise  ratio  of  2L 
over  that  obtained  with  a single  impulse. 

Golay  codes 

Consider  a digital  signal  composed  of  a series  of  sample 
values,  a(k)  and  b(k) . Golay  codes  can  be  computed  from  an 
initial  pair,  a1#  bx : 

ai=  {+1  +l} 

bx=  {+1  -1} 

Extensions  of  these  two  sequences  can  be  generated  by 
appending  h1  to  a1  to  obtain  a2  and  appending  -bx  to  a±  to 
obtain  b2 . Thus,  the  second  Golay  sequence  is 

a2=  {+1  +1  +1  -1}  and  b2=  {+1  +1  -1  +1} 

This  procedure  may  be  repeated  recurrently  to  produce  code 
pairs  of  any  desired  length,  L=2N.  The  critical  property  of 
this  pair  of  codes  is  that  the  autocorrelation  of  each  code  is 
exactly  zero  everywhere  except  at  the  origin.  That  is,  let 
C (n)  be  the  sum  of  the  autocorrelation  function  of  the  two 


sequences 


19 


C(n)  = laN(j)  aN{n+j)  +bN(j)  bN(n+j)]  for-L<n<+L 
j 

Eq.  2-1 

where  a(i)  or  b(i)  is  the  ith  value  of  the  aN  or  bN  Golay  code 
and  n is  the  lag  value.  C(n)  is  exactly  zero  for  all  values 
of  n other  than  n=0,  where  the  value  is  2L.  Autocorrelation 
in  the  time  domain  can  be  replaced,  of  course,  by 
multiplication  of  complex  conjugates  in  the  frequency  domain. 
That  is, 

fft(aN)  fft*  {aN)  +fft(bN)  fft*  (bN)  =2L  Eq.  2-2 

where  fft  stands  for  the  fast  Fourier  transform,  and  fft*(aN) 
and  fft*(bN)  are  the  complex  conjugates  of  fft  (aN)  and  f ft  (bN)  . 

This  property  of  Golay  codes  can  be  employed  to  obtain 
the  system  description.  We  first  apply  aN  to  the  system  and 
obtain  fft(aN)H(w)  as  the  output,  where  H(w)  is  the  system 
frequency  response,  the  Fourier  transform  of  the  impulse 
response.  Next,  we  apply  bN  to  the  system  and  obtain 
fft(bN)H(w)  as  the  output.  Finally,  by  multiplying  each 
output  with  the  conjugate  of  the  appropriate  Golay  code,  we 
determine  the  system  frequency  response. 

[fft(aN)  H{  w)  ] fft*  ( aN ) + [ fft  ( bN ) ff(o>)  ] fft*  ( bN ) 

= [fft(aN)  fft*  (aN)  +fft(bN)  fft*  (bN)]H(co)  Eq.  2-3 

= 2 LH  ( to ) 

Signal -to-noise  ratio 

The  power  of  any  signal  is  equal  to  the  autocorrelation 
function  of  the  signal  evaluated  at  zero  lag.  For  an  impulse 
of  unit  height  the  signal  power  is  unity.  For  the  Golay  codes 
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of  length  L the  signal  power  is  2L,  as  we  have  discussed 
previously.  Recall  that  the  autocorrelation  function  for  a 
single  pulse  is  nonzero  only  at  zero  lag.  Thus,  the  power 
spectrum  of  a single  impulse  and  the  power  spectrum  of  Golay 
codes  are  identical  in  shape  but  differ  in  level  by  101og(2L) 
dB . Therefore,  for  any  background  noise,  the  difference  in 
signal-to-noise  ratio  at  any  frequency  is  also  10log(2L) . In 
our  measurement,  a pair  of  Golay  codes,  with  L = 512,  were 
used  to  measure  the  frequency  response  of  the  system. 
Compared  with  a single  impulse,  the  improvement  of  signal-to- 
noise  ratio  is  over  30  dB . 

Measurements  and  Discussion 

Frequency  response  of  the  physical  model  with  the  two  driving 
sources 

The  frequency  response  in  the  physical  model  with  both 
the  long  lossy  tube  and  the  hard-wall  material  as  the  driving 
source  was  recorded  by  an  electret  microphone.  The  microphone 
was  calibrated  against  a standard  Briiel  & Kjaer  microphone 
(model  number  4133)  . The  frequency  response  of  the  system  was 
measured  using  Golay  codes. 

Frequency  response  of  the  sound  sources.  To  estimate  the 
frequency  response  of  the  sound  sources,  the  long  lossy  tube 
(2  meter)  was  connected  to  an  even  longer  (5  meter)  plastic 
tube,  which  makes  an  essentially  reflectionless  termination, 
because  the  attenuation  produced  by  the  5 -meter  tube  was 
nearly  120  dB . The  microphone,  located  at  the  end  of  the  long 
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lossy  tube  (2  meter) , recorded  the  frequency  response  shown  by 
the  curve  in  the  top  panel  of  Fig.  2-3.  The  frequency 
response  of  the  long  lossy  tube  source  includes  the  frequency 
response  of  the  amplifier,  the  tweeter,  the  horn,  as  well  as 
the  long  lossy  tube. 

The  frequency  response  of  the  hard-wall  source  was 
measured  in  the  same  way  except  that  the  hard-wall  material 
was  inserted  just  before  the  microphone  (see  Fig.  2-2a) . The 
curve  in  the  bottom  panel  of  Fig.  2-3  shows  the  frequency 
response  measured  with  the  hard-wall  material  in  place.  It 
includes  the  frequency  response  of  the  long  lossy  tube  source 
plus  the  frequency  response  of  the  hard-wall  material . 

Frequency  response  at  the  position  of  the  eardrum  with  a 
long  lossy  tube  as  the  driving  source.  The  9.7  cm  physical 
model  was  driven  by  the  long  lossy  tube,  and  the  microphone 
was  located  at  the  position  of  the  eardrum- -the  right  most  end 
of  the  physical  model.  Subtracting  the  frequency  response  of 
the  long  lossy  tube  (top  panel  of  Fig.  2-3)  from  the  frequency 
response  measured  at  this  position,  the  frequency  response  of 
the  model  for  a fixed  input  level  is  shown  by  the  curve  in 
panel  A of  Fig.  2-4. 


RELATIVE  RESPONSE  FOR  FIXED  INPUT  LEVEL  (dB) 
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Figure  2-3.  The  top  panel  is  the  frequency  response  of  the 
long  lossy  tube  source.  The  bottom  panel  is  the  frequency 
response  of  the  hard-wall  source. 
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Figure  2-4.  The  curve  in  panel  A represents  the  frequency 
response  at  the  position  of  the  eardrum  with  fixed  input  level 
when  the  sound  source  is  the  long  lossy  tube.  Panel  B the 
frequency  response  of  the  physical  model  with  the  hard-wall 
material  as  the  driving  source.  The  solid  and  dotted  lines 
are  the  frequency  responses  recorded  at  the  positions  of  the 
entrance  and  the  eardrum,  respectively.  Panel  C shows  the 
comparison  of  the  frequency  response  of  a real  ear  and  the 
physical  model.  The  solid  line  is  the  frequency  response  of 
the  tube-earmold-canal  system  recorded  at  the  entrance.  The 
dotted  line  is  the  frequency  response  of  the  physical  model 
recorded  at  the  entrance. 
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At  the  end  of  the  physical  model  driven  by  a long  lossy  tube, 
there  is  very  little  change  in  sound  pressure  as  a function  of 
frequency  for  a fixed  level  of  input.  The  amplitude  is 
approximately  6 dB  because  of  pressure  doubling  at  the  hard 
termination.  If  the  physical  model  is  a valid  approximation 
of  the  tube-earmold-canal  system,  then  it  can  be  concluded 
that  a long  lossy  tube  delivering  a fixed  input  level  produces 
essentially  no  change  in  pressure  as  a function  of  frequency 
at  the  eardrum. 

Frequency  response  at  the  positions  of  the  entrance  and 
the  eardrum  with  a hard-wall  driving  source.  The  frequency 
response  of  the  physical  model  was  measured  at  two  positions 
when  the  driving  source  was  the  hard-wall  material.  Again  the 
frequency  response  of  the  hard-wall  source  (bottom  panel  of 
Fig.  2-3)  was  subtracted  from  the  frequency  responses  measured 
at  the  positions  of  the  entrance  and  the  eardrum.  The  solid 
and  dotted  lines  in  panel  B of  Fig.  2-4  show  the  frequency 
responses  recorded  at  the  entrance  and  the  eardrum, 
respectively.  Note  that  the  resonance  and  anti-resonance 
frequencies  are  the  same  for  both  locations.  The  levels  of 
the  maximum  sound  pressure  are  the  same  for  both  locations, 
but  the  levels  of  the  minimum  sound  pressure  are  different. 
As  will  be  discussed  later,  the  ratio  between  the  maximum  and 
the  minimum  sound-pressure  level  at  the  positions  of  the 
entrance  and  the  eardrum  is  approximately  2 to  1 on  a 
logarithmic  scale . 
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It  is  important  to  note  that  the  patterns  seen  in  panel 
B of  Fig.  2-4  resemble  standing  wave  patterns,  but  are 
actually  variations  in  pressure  as  a function  of  frequency. 
The  standing  wave  pattern  describes  the  variation  of  sound 
pressure  as  a function  of  location  for  a fixed  frequency. 
Both  arise  because  of  resonance  in  an  acoustic  tube.  The 
focus  of  this  paper  is  on  the  variation  of  sound  pressure  as 
a function  of  frequency  for  a fixed  location.  Presumably  the 
variation  of  this  pressure  at  the  eardrum  influences  the 
changes  in  absolute  thresholds  measured  as  a function  of 
frequency . 

An  important  parameter  which  characterizes  the 
reflectance  of  an  acoustic  tube  is  called  standing  wave  ratio 
(SWR)  . It  is  defined  as  the  difference  between  the  maximum 
and  minimum  sound-pressure  levels  along  the  tube  for  a fixed 
input  frequency.  In  this  paper,  another  similar  parameter, 
the  maximum  to  minimum  sound  pressure  ratio  (MMSPR) , is 
defined  as  the  difference  between  the  maximum  and  minimum 
sound-pressure  levels  as  a function  of  frequency  at  a fixed 
location . 

Frequency  response  of  a real  ear  with  the  hard-wall  material 
as  the  driving  source 

To  demonstrate  that  the  physical  model  is  a valid 
approximation  of  the  real  ear,  the  frequency  response  was 
measured  with  a microphone  located  at  the  entrance  to  the 
tube-earmold-canal  system.  The  solid  line  in  panel  C of  Fig. 
2-4  shows  the  measured  pressure  for  one  ear.  The  frequency 


26 


response  of  the  physical  model  when  recorded  from  the  same 
position  is  displayed  as  the  dotted  line.  These  curves 
demonstrate  how  sound  pressure  changes  as  a function  of 
frequency  for  fixed  input  levels.  Note  that  the  patterns  are 
similar  for  the  real  ear  and  the  physical  model.  MMSPRs  are 
smaller  in  the  real  ear  data  than  in  the  physical  model.  This 
is  because  the  eardrum  impedance  is  much  lower  than  the 
impedance  of  the  hard  aluminum  employed  in  the  physical  model. 

The  objective  of  the  above  measurements  was  to  estimate 
the  sound  pressure  at  the  eardrum.  The  approach  was  to 
measure  the  pressure  at  the  entrance  to  the  ear  canal,  as 
shown  by  the  solid  line  in  panel  C of  Fig.  2-4,  and  use  a 
single  physical  model  to  infer  the  pressure  near  the  eardrum. 
The  relationship  between  the  pressure  measured  at  the 
positions  of  the  entrance  and  the  eardrum  can  be  observed  in 
panel  B of  Fig.  2-4.  Next,  equations  that  specify  the  sound 
pressure  level  at  various  positions  in  the  physical  model  will 
be  derived. 

Sound  pressure  variation  with  frequency  at  various  locations 
in  a closed  tube 

General  solution  of  the  sound  pressure  in  the  physical 
model . The  physical  model  is  shown  in  Fig.  2-2b.  The  uniform 
cylindrical  tube  is  closed  at  one  end  by  a wall  with  an 
acoustic  impedance  Ze,  and  fitted  with  a rigid,  frictionless 
sound  source  at  the  other  end.  In  the  measurement,  the  source 
is  a hard-wall  material,  and  can  be  regarded  as  a constant 
velocity  source.  It  moves  according  to 
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up=u0c os  (oat)  Eq.  2-4 

where  up  is  the  velocity  of  the  piston  and  u0  is  the  velocity 
amplitude . 

For  frequency  f < 0.586(c0/d),  where  c0  is  the  velocity 
of  the  sound  and  d is  the  diameter  of  the  tube,  only  the 
plane-wave  mode  can  propagate  in  the  tube  (Stevens  et  al . , 
1987)  . In  the  measurement,  d = 0.8  cm  and  c0  = 34300  cm/s,  so 
for  f < 25000  Hz,  the  sound  waves  in  the  tube  can  be 
considered  as  purely  longitudinal  plane  waves.  Therefore,  the 
steady  sound  field  in  the  tube  satisfies  the  one-dimensional 
plane  wave  equation 

d2u  _ 1 d2u 

d2x  c02  d2t  Eq>  2~5 

where  u is  the  velocity  of  the  medium. 

The  general  solution  to  Eq.  2-5  takes  the  form 

u{x,  t)  =A0eUkx-u,t) +B0e-Ukx+u,t)  Eq.  2-6 

where  u(x,t)  is  the  velocity  of  the  medium  at  a specific 
position  x and  moment  t,  k is  the  wave  number  equal  to  27rf/c0, 
and  f is  the  frequency  of  the  sound.  A0  and  B0  are  two 
constants  determined  by  boundary  conditions.  The  first  term 
of  Eq.  2-6  represents  the  wave  traveling  down  the  tube,  and 
the  second  term  represents  the  wave  reflected  back. 

As  the  acoustic  impedance  at  the  end  of  the  tube  is  Ze, 
the  first  boundary  condition  is 
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p(L,  t) 
u{L,  t) 


Eq.  2-7 


The  second  boundary  condition  is 


u( 0,  t)  =u0c os  (tot) 


Eq.  2-8 


in  which  u0cosojt  represents  the  velocity  of  the  piston. 

Using  these  two  boundary  conditions,  A0  and  B0  in  Eq.  2-6 
can  be  solved.  The  sound  pressure  is 


p {x,  t ) =fcos  (oj  t) 


( PqCq  + Ze)  e~lk{L~x)  - ( p0c0-Ze) 
(Poco  + Ze)  e'ikL+  (Poco"^e) 


q ik{L-x) 

elkL 


Eq.  2-9 


where  \p  is  the  constant  p0c0u0.  Here  p0  is  the  density  of  the 
medium . 


If  a plane  wave  traveling  in  a fluid  has  an  amplitude  (3 
at  position  x0,  then  its  amplitude  at  position  x will  be 


Pe 


-a  (x-x0) 


Eq.  2-10 


The  quantity  a is  known  as  the  attenuation  coefficient,  and 
represents  the  loss  in  the  tube  when  sound  travels  from  x0  to 
x due  to  viscosity.  This  loss,  omitted  in  the  above 
discussion,  may  be  incorporated  in  the  solution  by  allowing  k 
in  Eq.  2-6  to  be  a complex  number.  Setting  k = k + ia,  the 
new  general  solution  is 
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u(x,  t)  =A0e-axei{kx~u>t) +B0eaxe-iikx+wt)  Eq.  2-11 

The  solution  for  the  pressure  waveform  can  be  found  by  simply 
setting  k = k+ia  in  Eq.  2-9.  The  sound  pressure  at  the 
position  of  the  entrance  (x  = 0 ) is  then 


P (0  , t)  =lJfCOS  (co  t) 


(p0c0  + Ze)  eaL~ikL-  (p0c0-Ze)  e~aL+ikL 
(p 0c0+Ze)  eaL~ikL+  (p 0c0-Ze)  e~aL+lkL 


Eq.  2-12 


and  the  sound  pressure  at  the  position  of  the  eardrum  (x  = L) 
is 


p(L,  t ) = i|xcos  ( a)  t) 


2 Za 


(Poco  + Ze'>  e 


aL-ikL 


+ (PoC0~Ze'>  e 


-aL+ikL 


Eq.  2-13 


For  convenience,  in  the  following  discussion,  the 
termination  reflectance  Re  is  defined  as  (Ze  - p0c0)  / (Ze  + p0c0)  . 
The  magnitude  of  the  sound  pressure  at  the  entrance  and  the 
eardrum  are 


p(x=  0)  =201og|i|j 


Q<3.L-ikL+  g Q-aL+ikL 
QO.L-ikL_R  g-aL+ikL 


Eq.  2-14 


and 


p(x=L)  =201og|i|r 


1 +Rr 


gaL-ikL_j^  Q-aL+ikL 


Eq.  2-15 


MMS PR -Maximum  to  minimum  sound  pressure  ratio.  For  the 
current  investigation,  it  is  important  to  know  the  maximum  to 
minimum  sound  pressure  ratio  (MMSPR)  at  the  positions  of  the 
entrance  and  the  eardrum.  The  MMSPR  is  the  difference  between 
sound-pressure  levels  at  frequencies  of  resonance  and  anti- 
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resonance.  It  reflects  the  size  of  change  in  sound  pressure 
as  a function  of  frequency  at  a specific  position.  The 
resonance  and  anti-resonance  points  can  be  found  by  taking  the 
derivatives  of  Eqs . 2-14  and  2-15,  as  shown  in  Table  2-1. 

Table  2-1  shows  that  the  MMSPR  is  a function  of  the 
magnitude  of  reflectance  a , the  attenuation  coefficient  oi,  and 
the  length  of  the  tube  L.  It  is  independent  of  the  phase  of 
the  reflectance  0.  What  is  interesting  is  that  the  MMSPR  at 
the  entrance  and  the  eardrum  differ  exactly  by  a factor  of  2 
on  a logarithmic  scale. 

Estimation  of  g.  The  attenuation  coefficient  of  the 
plastic  tube,  a,  was  measured  as  follows.  A small  microphone 
was  placed  at  two  locations  in  the  plastic  tube  100  cm  apart, 
and  a transfer  function  was  measured  at  each  location.  The 
difference  between  the  transfer  functions  represents  the 
attenuation  of  sound  in  the  plastic  tube  across  frequency.  As 
might  be  expected,  the  attenuation  increases  with  frequency. 
Figure  2-5  shows  the  measured  attenuation  coefficients  (open 
circles)  as  calculated  from  Eq.  2-10.  As  the  attenuation 
coefficient  increases  almost  linearly  with  frequency,  a best 
fitting  solid  line  (slope  = 0.076/kHz,  interception  = 0.53)  is 
used  to  represent  the  attenuation  coefficient  across 
frequency . 


Table  2-1:  Resonance  and  anti-resonance  in  the  physical  model 
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FREQUENCY  (kHz) 


Figure  2-5.  Attenuation  coefficients  of  the  plastic  tube. 
The  open  circles  are  calculated  attenuation  coefficients  of 
the  plastic  tube  across  frequencies  from  2 to  16  kHz  (Eq.  2- 
10)  . The  solid  line  is  a linear  approximation  to  the  circles. 
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Values  of  a and  6.  Hudde  (1983)  reported  that  the 
estimates  of  the  eardrum  impedance  varied  widely  across 
different  investigations.  Individual  differences  in  ear  canal 
morphology  and  differences  among  measurement  techniques  may 
contribute  to  the  different  estimates  of  the  eardrum 
impedance . 

Two  different  definitions  of  eardrum  impedance  result  in 
two  different  methods  of  measuring  eardrum  impedance.  One  is 
to  choose  a reasonable  point  on  the  eardrum  and  measure  the 
sound  pressure  level  and  eardrum  velocity  at  that  point . For 
high  frequencies,  however,  as  the  sound  pressure  is  different 
across  the  eardrum  (Stinson  1985) , it  is  almost  impossible  to 
find  a "reasonable"  point  on  the  eardrum  to  measure  the 
eardrum  impedance.  A second  method  is  to  define  an  "eardrum 
reference  plane"  in  front  of  the  eardrum,  and  measure  the 
sound  pressure  and  velocity  at  the  plane.  The  latter 
definition  does  not  reproduce  the  mechanical  impedance  of  the 
eardrum.  It  does,  however,  reproduce  the  effective  acoustic 
termination  of  the  ear  canal,  which  is  appropriate  for  a 
physical  model  of  the  ear  canal . 

Using  the  second  method,  Hudde  (1983)  measured  the  sound 
pressure  at  various  locations  in  the  ear  canal  with  a probe 
microphone.  From  these  measurements,  he  deduced  the 
reflectance  of  the  eardrum  for  six  people.  These  results  are 
reproduced  as  the  solid  lines  in  Fig.  2-6.  The  dotted  curve 
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Figure  2-6.  Magnitude  and  phase  of  the  eardrum  reflectance. 
The  solid  lines  are  the  average  magnitude  and  phase  of  the 
eardrum  reflectance  for  six  individuals  from  Hudde  (1983) . 
The  dotted  line  in  the  top  panel  is  a third-order-polynomial 
approximation  to  the  magnitude.  The  dotted  line  in  the  bottom 
panel  is  a linear  approximation  to  the  phase. 
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in  the  top  panel  is  a third-order-polynomial  approximation  of 
the  magnitude  of  eardrum  reflectance  a.  As  the  phase  of  the 
reflectance  0 is  an  approximately  linear  function  of 
frequency,  the  straight  dotted  line  in  the  bottom  panel  is 
used  to  represent  0.  The  recent  study  by  Voss  and  Allen 
(1994)  on  acoustic  reflectance  in  the  ear  canal  showed  that 
there  were  large  across-subject  differences  in  the 
reflectance.  Their  average  result  agreed  reasonably  well  with 
the  data  of  Stinson  (1990)  and  Hudde  (1983)  . 

Calculations  of  the  sound  pressure  in  the  physical  model. 
Using  Hudde' s values  of  a and  0,  and  the  estimated  values  of 
oc  in  Fig.  2-5,  Eqs . 2-14  and  2-15  were  used  to  calculate  the 
sound  pressure  at  the  positions  of  the  entrance  and  the 
eardrum.  The  calculated  sound  pressures  are  shown  in  Fig.  2-7 
by  the  solid  (at  the  entrance)  and  dotted  (at  the  eardrum) 
lines,  respectively.  In  the  calculation,  the  length  of  the 
tube  L was  set  at  9.7  cm.  In  the  figure,  the  maximum  and 
minimum  sound  pressure  at  the  entrance  (solid  line)  shows 
little  change  across  frequency;  on  the  contrary,  the  maximum 
and  minimum  sound  pressure  at  the  eardrum  (dotted  line)  drifts 
downward  by  more  than  15  dB  over  the  frequency  range  of  5 to 
16  kHz. 

The  theoretical  derivation  above  demonstrates  that  the 
sound  pressure  at  the  entrance  and  at  the  eardrum  can  be 
predicted  by  Eqs.  2-14  and  2-15,  respectively.  Because  the 
value  of  ot  can  be  estimated  with  fair  accuracy,  there  are 
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FREQUENCY  (kHz) 


Figure  2-7.  Sound  pressure  at  the  positions  of  the  entrance 
and  the  eardrum  in  the  physical  model.  The  values  in  Fig.  2-5 
were  used  as  the  eardrum  (termination)  reflectance.  The  solid 
line  is  the  sound  pressure  at  the  entrance,  and  the  dotted 
line  is  the  sound  pressure  at  the  eardrum. 
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essentially  three  free  parameters  in  both  equations,  L,  a , and 

4>. 

Summary . The  tube-earmold-canal  system  can  be 
approximated  by  a physical  model  with  a uniform  cross- 
sectional  area.  The  frequency  response  of  the  physical  model 
with  both  the  long  lossy  tube  and  the  hard-wall  material  as 
the  driving  source  was  measured.  It  was  found  that  there  were 
no  changes  in  pressure  as  a function  of  frequency  at  the 
position  of  the  eardrum  when  the  ear  canal  was  driven  by  the 
long  lossy  tube.  There  was,  however,  appreciable  sound 
pressure  variation  as  a function  of  frequency  at  the  position 
of  the  eardrum  when  the  ear  canal  was  driven  by  the  hard-wall 
source.  The  equations  that  relate  the  sound  pressure  at 
various  locations  in  the  physical  model  were  derived.  Thus, 
if  the  parameters  can  be  correctly  estimated  from  the  measured 
sound  pressure  at  the  entrance  to  the  ear  canal,  the  sound 
pressure  at  the  eardrum  can  be  predicted  using  these  estimated 
parameters . 


Estimation  of  Pure-tone  Threshold 

In  the  following  two  experiments,  absolute  pure-tone 
thresholds  at  a variety  of  frequencies  were  measured  using 
both  the  hard-wall  and  the  long  lossy  tube  source.  As 
discussed  earlier,  the  different  sources  cause  different 
variations  in  sound  pressure  with  frequency  at  the  eardrum. 
The  objective  of  this  experiment  was  to  determine  whether 
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these  changes  in  pressure  produced  corresponding  changes  in 
absolute  threshold. 

Experiment  1 

Subi ects 

All  three  subjects  were  students  at  the  University  of 
Florida.  They  were  recruited  from  an  advertisement  in  the 
campus  newspaper,  and  were  paid  for  their  participation.  The 
subjects  were  hired  based  upon  the  following  criteria.  First, 
the  cross-sectional  area  of  the  subject's  ear  canal  was 
approximately  0.8  cm  in  diameter,  and  the  tragus  did  not 
prevent  easy  access  to  the  ear  canal  with  an  earmold.  Second, 
the  impulse  response  measured  at  the  entrance  to  the  ear  canal 
had  clear  resonance  and  ant i -resonance  features,  and  the 
frequency  spacing  was  similar  to  that  seen  in  a simple  uniform 
tube.  Third,  each  subject's  audiogram  when  measured  with  a 
Bekesy  audiometer  was  within  15  dB  of  normal  hearing  (ISO 
standard)  in  the  range  from  125  to  8000  Hz. 

Each  subject's  ear  canal  was  cleaned  in  the  university 
infirmary,  and  the  subject's  earmold  was  placed  in  the  ear 
canal  before  any  threshold  measurements  were  taken. 

Stimulus  and  procedure 

The  listeners  were  tested  in  a double-walled  booth.  The 
same  procedures  were  used  for  both  the  long  lossy  tube  and  the 
hard-wall  source.  First,  the  earmold  and  short  plastic  tube 
were  attached  and  positioned  in  the  pinna.  This  system  was 
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also  held  in  position  by  headgear  to  minimize  changes  in 
position  between  the  listener's  head  and  the  earmold.  A 
maximum- likelihood  procedure  was  used  to  determine  the 
absolute  threshold  in  a two-alternative,  forced-choice  task. 
Frequencies  between  5000  and  16,000  Hz  were  tested  in  200-Hz 
steps.  The  signal  duration  was  250  ms.  The  threshold  was 
estimated  as  the  signal  level  corresponding  to  the  94%  point 
on  the  most  likely  psychometric  function  (Green,  1990)  . The 
threshold  estimate  for  each  frequency  was  taken  as  the  average 
of  a session  of  seven  blocks  of  30  trials. 

To  insure  that  the  position  of  the  source  did  not  change 
during  the  course  of  the  experiment  and  to  insure  that  the 
source  was  placed  at  approximately  the  same  distance  from  the 
eardrum,  frequent  measurements  were  made  of  the  impulse 
response  of  the  system.  The  impulse  response  (with  a 
frequency  resolution  of  about  100  Hz)  was  recorded  with  a 
microphone  mounted  at  the  entrance  to  the  tube-earmold-canal 
system  at  the  beginning  and  end  of  each  session.  If  the  RMS 
difference  between  the  two  impulse  responses  was  greater  than 
1.5  dB  in  the  frequency  range  from  5000  to  16,000  Hz,  the  data 
for  that  session  were  discarded.  The  earmold  was  removed 
after  one  session,  and  the  subject  took  a ten-minute  break. 
At  the  beginning  of  a second  session,  the  position  of  the 
tube-earmold-canal  system  was  adjusted  so  that  the  RMS  of  the 
impulse  response  was  within  1.5  dB  of  that  recorded  at  the  end 
of  the  previous  session. 
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Figure  2-8  demonstrates  the  repeatability  of  frequency 
responses  at  the  position  of  the  entrance  across  different 
sessions  of  threshold  measurements.  This  figure  shows  the 
average  frequency  response  measured  at  the  beginning  of  five 
different  sessions  for  one  subject.  The  error  bars  represent 
standard  deviations  at  different  frequencies  across  the  five 
sessions.  The  average  standard  deviation  is  1.1  dB . 

Results  and  discussion 

Threshold  results.  Pure- tone  thresholds  (dB  sound 
pressure  level)  for  three  subjects  are  displayed  in  Fig.  2-8 
for  the  frequency  range  from  5000  to  16,000  Hz.  The  top 
curves  are  measured  thresholds  with  the  long  lossy  tube  as  the 
driving  source.  The  error  bars  are  the  standard  errors  across 
seven  sessions.  The  thresholds  represent  pressure  levels 
measured  at  the  entrance  of  a reflectionless  termination 
connected  to  the  long  lossy  tube  source.  The  pressure  at  the 
eardrum  should  be  6 dB  higher  than  the  measured  pressure  at 
the  entrance  (see  panel  A in  Fig.  2-4) . The  thresholds  with 
the  long  lossy  tube  as  the  driving  source  reflect  the  change 
in  sensitivity  of  the  listener  as  a function  of  frequency. 

The  middle  set  of  curves  were  obtained  with  the  hard-wall 
source.  The  thresholds  represent  pressure  levels  measured  at 
the  entrance  of  a reflectionless  termination  connected  to  the 
hard-wall  source.  These  pressures  reflect  changes  in  the 
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Figure  2-8.  The  average  frequency  response  at  the  entrance 
recorded  in  five  different  sessions  of  threshold  measurements. 
The  error  bars  are  the  standard  deviations  across  the 


sessions . 
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Figure  2-9.  Threshold  results  for  three  subjects.  The  top 
panels  are  thresholds  in  dB  SPL  with  the  long  lossy  tube  as 
the  driving  source.  The  middle  panels  are  thresholds  in  dB 
SPL  with  the  hard-wall  material  as  the  driving  source.  The 
bottom  panels  are  the  difference  between  the  middle  and  top 
panels . 
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sensitivity  of  the  listener  as  a function  of  frequency  as  well 
as  the  changes  in  sound-pressure  level  at  the  eardrum  caused 
by  the  standing  waves  created  by  the  hard-wall  source.  The 
effect  of  the  standing  waves  can  be  determined  by  subtracting 
the  curves  of  the  top  panels  from  those  of  the  middle  panels. 
This  difference  is  shown  in  the  bottom  panels  of  the  figure. 
It  is  clear  that  the  changes  in  sound  pressure  at  the  eardrum 
produce  the  changes  in  threshold  as  a function  of  frequency. 
The  differences  between  the  peaks  and  the  valleys  are  about  8- 
15  dB. 

In  the  next  section,  the  sound  pressure  variation  with 
frequency  at  the  eardrum  will  be  derived,  and  this  function 
will  then  be  compared  with  resulting  threshold  variations 
shown  in  the  bottom  panels  of  Fig.  2-9.  Before  doing  that,  we 
compare  the  present  threshold  measurements  with  those  obtained 
by  previous  investigators. 

The  estimated  pressure  at  the  eardrum  corresponding  to 
the  thresholds  measured  with  the  long  lossy  tube  are  plotted 
as  filled  circles  in  Fig.  2-10.  They  represent  the  average  of 
the  three  subjects.  Threshold  estimates  from  several  previous 
studies  are  also  plotted  in  the  same  figure.  Large 
variability  was  expected  across  the  studies  as  different 
sound-delivery  systems  and  calibration  procedures  were  used  by 
different  investigators.  Four  of  the  six  studies  show  little 
increase  (=  10  dB)  in  absolute  threshold  as  the  signal 
frequency  changes  from  8 to  16  kHz.  Two  of  the  six  studies 
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Figure  2-10.  Comparison  of  thresholds  estimated  in  this  study 
(solid  circles)  and  those  of  previous  studies. 
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show  a much  large  change  (=  30  dB) . The  thresholds  obtained 
in  the  present  study  are  generally  similar  (±15  dB)  to 
thresholds  measured  in  previous  studies. 

Model  fitting.  If  the  ear  canal  is  treated  as  a tube 
with  a uniform  cross-sectional  area,  then  the  sound  pressure 
at  the  eardrum  can  be  derived  from  the  sound  pressure  measured 
at  the  entrance  (Eqs.  2-14  and  2-15).  The  relationship 
depends  on  three  free  parameters,  the  eardrum  reflectance  Re, 
a complex  number  expressed  by  a magnitude  ,cr,  and  a phase  ,0, 
and  the  distance  from  the  sound  source  to  the  eardrum  L.  The 
general  approach  is  to  adjust  Re  and  L in  Eq.  2-14,  so  as  to 
minimize  the  rms  deviation  between  the  measured  and  predicted 
values  of  pressure  levels  at  the  entrance  to  the  ear  canal. 
Once  the  best  fitting  values  of  Re  and  L are  determined,  Eq. 
2-15  can  be  used  to  predict  the  expected  response  at  the 
eardrum,  P (x=L) . 

As  L is  small  (=  10  cm)  , the  impact  of  a in  solving 
Eqs.  2-14  and  2-15  is  insignificant.  The  measured  value  of  a 
taken  from  the  plastic  tube  (Fig.  2-5)  is  used  as  the 
attenuation  coefficient  of  the  tube-earmold-canal  system. 

The  key  issue  is  to  determine  the  eardrum  reflectance  Re . 
The  average  values  obtained  by  Hudde  (1983)  cannot  be  used, 
because,  as  he  noted,  there  are  large  individual  differences 
among  subjects.  His  average  data  are  representative  of  some 
average  eardrum  reflectance,  but  are  not  useful  in  trying  to 
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predict  the  individual  data  of  the  seven  listeners  in  this 
study . 

(1)  Estimation  of  a 

Table  2-1  shows  that  given  a value  of  L and  ot,  the 
magnitude  of  the  eardrum  reflectance  a is  expressed  by 

MMSPR 

a_c2ttL  ( 10  40  -1)  Eq.  2-16 

MMSPR 

(10  40  +1) 


Thus,  from  measurements  of  the  pressure  at  the  entrance  to  the 
ear  canal,  the  difference  between  successive  maximum  and 
minimum  values  can  be  used  to  determine  a different  value  for 
a using  Eq.  2-16.  A set  of  such  values  measured  at  different 
frequencies  can  be  approximated  by  a third  order  polynomial 
function,  so  that  the  variation  of  the  magnitude  of  the 
reflectance  changes  smoothly  with  frequency  cr  ( f)  . 

(2)  Estimation  of  0 and  L 

Estimation  of  the  phase  of  the  reflectance  is 
considerably  difficult,  because  of  a potential  coupling 
between  the  values  of  0 and  L.  To  illustrate  this  problem, 
assume  that  0 changes  linearly  with  frequency,  which  is 
approximately  true,  as  shown  by  Hudde  (1983)  . Table  2-1 
shows  that  resonance  frequency  of  the  tube  equals  [ (2n7T- 
0)  C]  / (4ttLi)  , {n  = 0,  1,  2,  ...}.  Thus,  the  frequency  spacing 
between  adjacent  resonance  frequencies  is 


A f= 


c 

2 L 


<i-M 

271 


) 


Eq.  2-17 


where  Af  is  the  frequency  spacing  between  adjacent  resonances, 
and  A 0 is  a measure  of  the  slope  of  the  phase  function  as  a 
function  of  frequency.  It  is  convenient  to  define  A 0 as  the 
difference  in  phase  measured  between  two  successive 
resonances.  If  0 is  a linear  function  of  frequency,  then  A 0 
is  constant  and  independent  of  frequency.  Therefore,  Eq.  2-17 
can  be  expressed  as 

Af= — — Eq.  2-18 
2 L* 

where  L*,  equal  to  L/  (l-A0/27r)  , is  the  "effective  length"  of 
the  tube.  For  a tube  of  length  L*,  the  frequency  spacing 
between  adjacent  resonances  is  the  same  as  Eq.  2-18  if  the 
reflectance  phase  of  the  tube  0*,  is  constant  and  independent 
of  frequency  (A 0*  = 0)  . 

To  see  the  effects  of  this  ambiguity,  the  pressure  at  the 
eardrum  is  calculated,  given  two  identical  measures  of 
pressure  at  the  entrance,  but  based  on  two  different 
assumptions  about  the  phase  of  the  reflectance.  Following 
Hudde  (1983),  assume  a listener  has  the  same  eardrum 
reflectance  as  the  average  values  obtained  from  his  average 
data  (dotted  lines  in  Fig.  2-6)  . The  change  in  0 over  two 
consecutive  resonances,  A 0,  has  a value  of  -0.56  radians.  If 
L is  equal  to  9.7  cm,  the  sound  pressure  measured  at  the 
entrance  is  given  by  the  solid  line  in  the  top  panel  of  Fig. 
2-11.  The  effective  length  L*  is  8.9  cm  (see  Eq.  2-18)  . The 
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Figure  2-11.  The  solid  line  in  the  top  panel  shows  the 
pressure  at  the  entrance  using  Hudde's  values  of  0 and  L = 9.7 
cm  (A 0 = -0.56) . The  dotted  line  shows  the  pressure  at  the 
entrance  with  0*  being  a constant  value  of  1.52  radians  and  L* 
= 8.9  ( A0*  = 0)  . The  curves  in  the  bottom  panel  show  the 
sound  pressure  at  the  eardrum  using  the  same  two  sets  of  0 and 
L. 
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dotted  line  in  the  top  panel  of  Fig.  2-11  shows  the  calculated 
sound  pressure  at  the  entrance  if  the  tube  has  an  effective 
length  of  8.9  cm  and  a constant  reflectance  phase  0 * = 1.52 
radians  (A0*  = 0)  . The  constant  reflectance  phase  will  be 
called  the  "equivalent  phase"  of  the  eardrum  reflectance.  The 
dotted  and  solid  curves  are  nearly  identical.  These  curves 
demonstrate  that  different  combinations  of  L and  0 can  create 
essentially  the  same  sound  pressure  at  the  entrance.  But  what 
is  the  pressure  at  the  eardrum  suggested  by  these  two  sets  of 
parameters? 

Two  sets  of  curves  can  be  calculated  for  the  pressure  at 
the  eardrum,  corresponding  to  these  two  assumptions  about  the 
character  of  0.  The  solid  line  in  the  bottom  panel  of  Fig.  2- 
11  is  the  sound  pressure  at  the  eardrum  calculated  from  Eq.  2- 
15,  using  Hudde's  values  as  the  eardrum  reflectance,  L = 9.7 
cm,  A0  = -0.56.  The  dotted  line  in  the  bottom  panel  of  Fig. 
2-11  is  calculated  from  Eq.  2-15,  using  L*  = 8.9  cm,  A0*  = 0 
(0*  =1.52  radians) . For  frequencies  lower  than  11  kHz,  there 
are  only  small  discrepancies  between  the  two  curves.  For 
frequencies  higher  than  13  kHz,  however,  the  discrepancy  can 
be  as  large  as  15  dB . 

Estimating  the  parameters  L and  0 from  the  pressure 
measured  at  the  entrance  to  the  ear  canal  cannot  provide 
unique  predictions  of  the  pressure  at  the  eardrum.  Serious 
discrepancies  can  arise  at  frequencies  higher  than  13  kHz. 
Note  that,  despite  this  failure  to  predict  the  same  sound- 
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pressure  levels,  the  frequencies  at  which  maximum  and  minimum 
pressure  occur  remain  nearly  the  same. 

(3)  Prediction  of  sound  pressure  at  the  eardrum 

This  ambiguity  as  to  the  phase  of  the  eardrum  reflectance 
makes  it  impossible  to  estimate  accurately  the  pressure  at  the 
eardrum  above  13  kHz.  We  dealt  with  this  problem  by  assuming 
that  cp  is  constant  and  independent  of  frequency,  A0  = 0 . The 
parameters  L*,  <f>* , and  o are  estimated  from  the  pressure 
measured  at  the  entrance  to  the  ear  canal.  The  rms  deviation 
between  the  predicted  (Eq.  2-14)  and  measured  sound  pressure 
at  the  entrance  is  minimized  by  adjusting  L*  and  <f>* . The  new 
value  of  L*  is  used  in  Eq.  2-16  to  produce  a new  cubic 
polynomial  function  a (f)  . The  procedure  is  repeated  until  the 
rms  error  reaches  a minimum  value.  The  typical  rms  value  was 
about  4 dB . 

Figure  2-12  shows  the  results  of  this  estimation 
procedure  applied  to  all  three  listeners.  The  predicted 
pressure  (Eq.  2-14)  at  the  entrance  to  the  ear  canal  is  shown 
as  the  dotted  lines,  whereas  the  open  circles  are  the  measured 
pressures.  The  rms  error  is  about  4 dB,  and  the  exact  value 
is  given  in  the  figure  along  with  the  best  fitting  values  for 
L*  and  . 

Figure  2-13  shows  the  magnitude  of  the  reflectance 
values,  a,  calculated  for  each  listener  and  the  polynomial 
fit  to  the  individual  data.  The  average  value  of  a for  the 
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Figure  2-12.  The  open  circles  are  the  measured  sound 
pressures  at  the  entrance  to  the  ear  canal.  The  dotted  lines 
are  the  optimum  fits  to  the  circles  obtained  from  Eq.  2-14. 
The  values  of  L*,  the  effective  length  of  the  physical  model, 
and  </>* , the  equivalent  phase  of  eardrum  reflectance,  for  the 
best  fits  are  shown  for  each  listener.  The  rms  deviation 
between  the  circles  and  the  dotted  line  is  also  shown  for  each 
individual . 


52 


5 7 9 11  13  15  5 7 9 11  13  15 

FREQUENCY  (kHz) 


Figure  2-13.  Estimates  of  a.  The  open  circles  are  the 
magnitude  of  eardrum  reflectance  between  successive  maximum 
and  minimum  pressure  points  calculated  from  Eq.  2-16.  The 
solid  curves  in  panels  A,  B,  and  C are  the  third-order- 
polynomial  approximations  of  a for  three  listeners.  The  solid 
curve  labeled  as  'mean'  is  the  average  a of  the  three 
listeners.  The  results  obtained  by  Stinson  (1990)  and  Hudde 
(1983)  are  also  displayed. 
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three  subjects  is  plotted  as  the  solid  line  labelled  as  'mean' 
in  the  figure.  The  estimates  of  the  magnitude  of  the  eardrum 
reflectance  in  this  study  are  about  20%  smaller  than  the 
results  obtained  by  Stinson  (1990)  and  Hudde  (1983) . 

With  the  estimated  values  of  L*,  0*,  and  a,  Eq.  2-15  is 
used  to  predict  the  value  of  pressure  measured  at  the  eardrum 
( P (x  = L) ) . Figure  2-14  shows  the  inverse  of  that  pressure 
along  with  the  difference  in  threshold  (replotted  from  the 
bottom  panels  of  Fig.  2-9) . A constant  has  been  added  to  the 
inverse  pressure  levels  to  best  fit  the  curve  to  the  observed 
data.  The  constant  reflects  the  absolute  sensitivity  of  the 
listener.  For  all  subjects,  the  data  appear  to  roughly  follow 
the  predicted  curve  below  13  kHz.  For  frequencies  beyond  13 
kHz,  possibly  because  of  the  poor  prediction  of  the  real  sound 
pressure  at  the  eardrum,  there  are  obvious  deviations  between 
the  predicted  curve  and  the  threshold  variation  for  subjects 
2 and  3.  Despite  these  differences,  the  number  of  maxima  and 
minima  in  the  predicted  curves  is  nearly  identical  to  that 
observed  in  the  threshold  variation  for  subjects  1 and  2. 
For  subject  3,  this  is  true  up  to  about  13  kHz.  Note  that, 
while  the  pressure  cannot  be  predicted  exactly  because  of  the 
phase  ambiguity,  the  number  of  maxima  and  minima  is  consistent 
with  estimated  values  and  Eq.  2-15. 

The  correlation  between  the  predicted  curve  and  the 
threshold  variation  is  0.53,  0.22,  0.1,  for  subjects  1,  2,  3, 
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Figure  2-14.  The  solid  curves  are  the  inverse  of  the 
predicted  sound  pressure  at  the  eardrum  calculated  from  Eq.  2- 
15.  The  solid  circles  are  the  threshold  differences  replotted 
from  the  bottom  panels  of  Fig.  2-9. 
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respectively.  The  small  correlation  coefficients  obtained  for 
subjects  2 and  3 are  largely  due  to  the  deviation  between  the 
predicted  curve  and  the  threshold  variation  above  13  kHz.  If 
only  the  frequency  range  below  13  kHz  is  considered,  the 
correlation  coefficients  increase  to  0.41  for  subject  2 and 
0.55  for  subject  3.  For  all  three  subjects,  in  the  frequency 
range  below  13  kHz,  the  correlation  coefficients  are 
significantly  different  from  zero  (p  < 0.01).  The  rms 
differences  between  the  predicted  curve  and  the  threshold  data 
below  13  kHz  are  3.7,  4.1  and  4.0  dB,  respectively.  For 
frequencies  higher  than  13  kHz,  the  discrepancies  are  3.3, 
6.8,  10.4  dB,  respectively. 

Because  of  the  coupling  of  0 and  L,  more  than  one  pair  of 
values  can  predict  exactly  the  same  variation  of  pressure  as 
a function  of  frequency  as  measured  at  the  entrance  to  the  ear 
canal . We  assumed  that  the  phase  0 was  constant  and 
independent  of  frequency  (A0=O) . This  assumption  was  used  to 
minimize  the  rms  difference  between  the  observed  variation  of 
pressure  as  a function  of  frequency  at  the  entrance  to  the  ear 
canal.  This  minimization  led  to  the  estimates  of  L,  a,  and  a 
used  in  Eq.  2-15  to  predict  the  variation  in  pressure  at  the 
eardrum.  These  predictions  correlate  only  modestly  with  the 
variation  in  threshold  that  we  actually  observed.  Probably 
some  other  assumptions,  in  particular  assuming  A0  * 0 , would 
produce  better  predictions.  One  possible  value  for  A0  is  - 
0.56,  the  average  slope  of  Hudde's  data  (Fig.  2-6).  Using 
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this  new  value  and  the  new  estimates  of  L,  a,  and  a,  Eq.  2-15 
was  once  more  used  to  predict  the  variation  of  pressure  as  a 
function  of  frequency  at  the  eardrum.  The  correlation  between 
these  pressure  variation  and  the  measured  thresholds  are  lower 
than  those  based  on  the  assumption  that  A0  = 0 . There  is  no 
point  in  generally  exploring  the  assumption  that  A0  ^ 0 . The 
set  of  values  are  unbounded,  and  every  possible  value  would 
produce  the  same  best  fit  to  the  pressure  variation  measured 
at  the  entrance  to  the  ear  canal.  In  fact,  the  variation 
would  be  exactly  the  same  as  that  produced  by  the  assumption 
of  A 0 = 0.  Thus,  given  only  those  data,  there  is  no 
principled  basis  to  choose  among  the  large  number  of 
equivalent  assumptions.  The  different  choices  can  only  be 
distinguished  on  the  basis  of  the  variation  in  pressure  as  a 
function  of  frequency  measured  at  the  eardrum,  which  is  what 
we  are  trying  to  predict.  Thus,  we  used  the  simplest 
assumption,  namely,  A0  = 0 . 

In  summary,  absolute  pure-tone  thresholds  were  measured 
from  5 to  16  kHz  with  and  without  the  influence  of  the 
standing  waves  in  the  ear  canal.  The  difference  in  these  two 
sets  of  thresholds  reflects  the  effect  of  changes  in  pressure 
at  the  eardrum  as  a function  of  frequency  created  by  standing 
waves.  Using  parameters  estimated  from  the  pressure  measured 
at  the  entrance  to  the  ear  canal,  the  pressure  at  the  eardrum 
was  calculated.  For  frequencies  below  13  kHz,  the  calculated 
sound  pressure  at  the  eardrum  predicted  the  changes  in 
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threshold 

with  an  accuracy  of 

about 

± 2 

dB . Above 

that 

frequency, 

the  error  increased 

to  as 

much 

as  ±5  dB . 

These 

large  errors  are  probably  caused  by  inaccurate  estimates  of 
the  phase  of  the  eardrum  reflectance  0. 

Experiment  2 

Subi ects 

Four  new  students  participated  in  the  experiment.  They 
were  screened  following  the  same  criteria  as  described  in 
Experiment  1 . 

Stimulus  and  procedure 

The  length  of  the  short  plastic  tube  which  connected  the 
earmold  and  the  hard-wall  material  was  reduced  from  7 cm  to 
approximately  3 cm.  Thus,  there  were  fewer  points  of  maxima 
and  minima  in  the  frequency  range  from  5000  to  16,000  Hz.  To 
find  the  largest  rms  difference  in  the  measured  threshold,  the 
test  frequencies  were  chosen  to  correspond  to  the  maximum  and 
minimum  pressure  points.  The  threshold,  with  and  without  the 
hard-wall  material,  for  four  new  subjects  was  tested  following 
the  same  procedure  described  in  experiment  1.  The  predicted 
sound  pressure  at  the  eardrum  for  each  individual  was  obtained 
using  the  model  fitting  procedure  described  previously. 
Results  and  discussion 

Differences  between  the  thresholds  with  the  long  lossy 
tube  source  and  the  hard-wall  source  are  shown  as  the 
asterisks  in  Fig.  2-15.  The  error  bars  represent  plus  or 
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Figure  2-15.  The  solid  curves  are  the  inverse  of  the 
predicted  sound  pressure  at  the  eardrum  for  four  new  subjects. 
The  asterisks  are  measured  threshold  differences;  error  bars 
indicate  the  standard  error  at  each  test  frequency.  The  rms 
difference  and  correlation  coefficient  between  the  asterisks 
and  the  predicted  data  at  the  test  frequencies  are  shown  for 
each  listener. 
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minus  one  standard  error.  The  inverse  of  the  predicted  sound 
pressure  at  the  eardrum  is  shown  by  the  solid  curve.  The 
correlation  coefficients  and  the  rms  differences  between  the 
measured  data  and  predicted  data  at  the  test  frequencies  are 
given  in  each  panel.  For  all  three  subjects,  in  the  frequency 
range  from  5000  to  16,000  Hz,  the  correlation  coefficients  are 
significantly  different  from  zero  at  the  0.05  level.  The  size 
of  the  correlation  indicates  that  changes  in  sound  pressure  at 
the  eardrum  do  cause  corresponding  changes  in  absolute 
threshold  for  the  test  frequency  region. 

Conclusions 

The  following  conclusions  can  be  drawn  from  the 
experiments . 

(1)  The  pressure  measurements  made  at  the  entrance  to 
the  ear  canal  with  a hard-wall  source  are  useful  in  predicting 
the  pressure  variation  at  the  eardrum.  For  a fixed  input 
level,  a good  estimate  of  the  change  in  pressure  as  a function 
of  frequency  at  the  eardrum  (termed  the  MMSPR)  is  one-half  the 
MMSPR  measured  at  the  entrance  to  the  ear  canal  on  a 
logarithmic  scale. 

(2)  The  calculated  changes  in  sound  pressure  at  the 
eardrum  can  predict  the  change  in  threshold  with  an  accuracy 
of  ±2  dB  below  13  kHz.  For  frequencies  above  13  kHz,  the 
error  can  be  as  large  as  ± 5 dB . 
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(3)  A long  lossy  tube  provides  a simple  means  of 
minimizing  the  changes  in  pressure  at  the  eardrum  created  by 
strong  resonances  and  anti-resonances  in  the  ear  canal. 

For  the  MMSPRs  measured  with  these  subjects,  changes  in 
pressure  at  the  eardrum  of  as  much  as  ± 5 dB  can  occur  if  the 
hard-wall  source  is  moved  ± 0.3  cm  at  13  kHz  (one-quarter 
wavelength) . Such  large  changes  in  pressure  do  not  occur  at 
low  frequencies  because  of  the  much  longer  wavelengths. 
Because  of  this  variation,  one  would  expect  the  variability  of 
threshold  estimates  to  be  much  larger  at  high  frequencies  than 
at  low  frequencies.  In  the  next  section,  we  will  investigate 
the  variability  of  intra- sub j ect  threshold  measurements  caused 
by  the  movement  of  transducers. 


Chapter  3 

RELIABILITY  OF  HIGH-FREQUENCY  THRESHOLD  MEASUREMENT 

Introduction 

In  this  chapter,  we  will  investigate  several  factors  that 
influence  the  reliability  of  high-frequency  threshold 
estimates.  As  discussed  by  Stelmachowicz  et  al . (1989),  the 

reliability  of  high-frequency  audiometry  is  determined  by  two 
variances.  The  first  is  called  threshold  variance.  For  any 
frequency,  two  threshold  estimates  obtained  within  the  same 
session  will  not  agree  perfectly,  even  if  the  measurement 
conditions  are  equivalent.  This  variance  is  determined 
largely  by  the  number  of  trials  tested,  the  slope  of  the 
underlying  psychometric  function,  the  stability  of  the 
observer,  and  the  psychophysical  procedure.  The  second  is 
called  fitting  variance . Due  to  the  existence  of  standing 
waves  in  the  ear  canal,  if  the  sound  source  is  removed  and 
replaced,  the  change  in  relative  position  between  the  sound 
source  and  the  eardrum  might  result  in  a change  of  sound 
pressure  in  the  ear  canal.  Thus  the  estimates  of  thresholds 
for  two  sound  source  positions  will  be  different.  This 
variance  is  determined  by  the  frequency  of  the  sound  source 
and  the  amount  by  which  the  sound  source  is  moved. 
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Three  experiments  were  designed  to  test  various  aspects 
of  this  problem.  In  the  first  experiment,  the  intra-subject 
threshold  variance  and  fitting  variance  were  measured  with  the 
hard-wall  source.  The  change  of  sound  pressure  in  the  ear 
canal  caused  by  the  movement  of  the  hard-wall  source  and  the 
resulting  variation  in  threshold  estimates  were  assessed.  In 
the  second  experiment,  we  estimated  the  change  of  sound 
pressure  in  the  ear  canal  caused  by  different  positioning  of 
headphones.  The  results  were  evaluated  in  terms  of  their 
effect  on  fitting  variance . This  was  done  for  two  types  of 
headphones  frequently  used  in  audiometric  testing- -one 
circumaural,  and  the  other  superaural . In  the  third 
experiment,  the  effect  of  the  slope  of  the  psychometric 
function  on  threshold  variance  was  investigated.  Previous 
experiments  suggested  that  the  slope  of  the  psychometric 
function  was  different  for  different  test  frequencies. 

Experiment  1 


Introduction 

In  Chapter  2,  it  was  estimated  that  sound  pressure  in  the 
ear  canal  changed  by  ± 5 dB  as  a function  of  frequency  because 
of  the  existence  of  standing  waves.  Similar  variation  of 
sound  pressure  can  be  produced  by  movement  in  the  hard-wall 
source  of  one -quarter  wavelength.  Will  the  movement  produce 
corresponding  variation  in  behavior  threshold  estimates.  In 
this  experiment,  the  variability  of  threshold  was  estimated  by 
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comparing  two  experimental  conditions.  In  the  first 
condition,  these  was  no  movement  in  the  position  of  the  hard- 
wall  source  across  all  test  sessions.  In  the  second 
condition,  the  position  of  the  hard-wall  source  changed 
between  each  test  session. 

Method 

Subi ects 

Four  subjects  with  normal  hearing  sensitivity  (250  to 
8000  Hz)  participated  in  the  experiment.  All  subjects  were 
students  at  the  University  of  Florida.  The  average  age  of  the 
subjects  was  21  years.  An  earmold  was  constructed  for  each 
subject  and  the  subjects  were  asked  to  have  their  external 
ears  cleaned  before  the  experiment.  Prior  to  data  collection, 
each  subject  completed  several  hours  of  practice  on  the  task. 
Subjects  were  paid  an  hourly  wage  and  received  a 20  percent 
bonus  upon  completion  of  the  study. 

Stimuli 

Low-  and  high-frequency  pure  tones  were  used  in  the  test 
to  compare  the  intra-subject  reliability  of  threshold 
estimates  for  both  frequencies.  To  find  the  largest  variation 
caused  by  the  movement  of  the  sound  source,  the  high  test 
frequencies  were  chosen  to  correspond  to  three  successive 
maximum  and  minimum  points  from  7 to  14  kHz  in  the  frequency 
response  of  the  ear  canal.  At  low  frequencies,  there  was 
little  variation  in  sound  pressure  as  a function  of  frequency. 
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and  only  1000  Hz  was  tested.  The  specific  test  frequencies 
for  each  subject  were  listed  in  the  second  column  of  Table  3- 
2 . 

Procedure 

The  pure-tone  signal  was  delivered  to  the  ear  canal  by 
the  hard-wall  source  as  described  in  Chapter  2.  For  each  test 
frequency,  the  pure-tone  threshold  was  measured  in  two 
conditions.  In  the  first  condition,  the  threshold  was 
measured  in  five  sessions  in  which  we  attempted  to  reposition 
the  sound  source  in  exactly  the  same  place  for  each  session. 
The  ability  to  achieve  this  goal  was  monitored  by  measuring 
the  frequency  response  of  the  ear  canal  before  and  after  each 
session.  For  any  data  reported  here,  the  rms  difference  in 
frequency  response  for  the  frequency  range  from  5 to  16  kHz 
was  always  less  than  1.5  dB . In  the  second  condition,  the 
threshold  was  measured  in  five  sessions  with  the  position  of 
the  hard-wall  material  changed  by  as  much  as  ± 0.5  cm  between 
two  sessions  of  measurement.  These  movements  produced  rms 
difference  in  the  frequency  response  greater  than  3 dB . 

For  each  session,  what  we  call  the  session  threshold  was 
an  average  based  upon  eight  blocks  of  30  trials.  A maximum- 
likelihood  procedure  was  used  to  determine  the  absolute 
threshold  in  a two-alternative,  forced-choice  task  ( 2AFC) . 
The  threshold  was  estimated  as  the  signal  level  corresponding 
to  the  94%  point  on  the  psychometric  function  (Green,  1990) . 
The  signal  duration  was  250  ms. 
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Results  and  Discussion 

The  standard  deviation  values  computed  across  the  five 
session  thresholds  are  displayed  in  Table  3-1  and  Table  3-2. 
Table  3-1  shows  the  result  for  a 1-kHz  signal  and  Table  3-2 
shows  the  result  for  high-frequency  signals.  The  standard 
deviation  values  shown  under  the  column  "without  movement" 
were  calculated  for  the  condition  in  which  the  distance 
between  the  hard-wall  source  and  the  eardrum  remained  the 
same.  The  standard  deviation  value  reflects  the  threshold 
variance  at  the  test  frequency.  The  standard  deviation  values 
shown  under  the  column  "with  movement"  were  obtained  for  the 
condition  in  which  the  distance  between  the  hard-wall  source 
and  the  eardrum  changed.  It  reflects  the  combination  effect 
of  both  threshold  variance  and  fitting  variance . 

An  ANOVA  is  used  to  test  the  null  hypothesis  that  the 
estimates  of  the  five  session  thresholds  are  the  same.  If  the 
null  hypothesis  can  be  rejected  at  the  0.05  level,  then  the 
computed  standard  deviation  is  marked  with  an  asterisk  (see 
Table  3-2)  . For  the  1-kHz  signal,  the  null  hypothesis  can  not 
be  rejected  for  all  the  points  tested.  The  threshold 
estimates  in  any  two  sessions  are  not  statistically  different 
from  each  other. 

For  high-frequency  signals,  when  the  sound  source  is  not 
moved,  the  null  hypothesis  can  be  rejected  for  only  three  of 
the  twelve  points  tested.  Perhaps  this  is  because  the 
criterion  of  1.5  dB  change  in  frequency  response  is  not 
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Table  3-1  Standard  deviation  of  1-kHz  signal  threshold  across 
five  test  sessions 


Without  Movement 

With  Movement 

Subj  ect 

STD  (dB) 

STD  (dB) 

SI 

0 . 65 

0 . 71 

S2 

0 . 79 

0 . 80 

S3 

0 . 88 

0 . 91 

S4 

0.43 

0.45 

AVG 

0 .69 

0 . 71 
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Table  3-2:  Standard  deviation  of  high-frequency  signal 

threshold  across  five  test  sessions 


Without  Movement 

With  Movement 

Subj  ect 

Freq (Hz ) 

STD  (dB) 

STD  (dB) 

SI 

9280 

0.36 

3 . 82* 

10740 

0 . 93* 

4.23* 

12110 

0 . 93* 

3.13* 

S2 

7420 

0 . 64 

2 . 70* 

9570 

0.34 

2 . 06* 

11040 

0.26 

3 . 75* 

S3 

9670 

0.41 

1.42* 

11130 

0 . 61 

2 . 83* 

12400 

0 . 83 

2 .69* 

S4 

10060 

0 . 80 

2 . 86* 

11820 

0 . 69 

3 . 23* 

13180 

1.09* 

1.41* 

AVG 

0 . 551 

2 . 85 

1.  Average  excluding  points  with  *. 
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sufficiently  stringent.  Subject  instability  is  probably 
another  reason.  These  three  points  were  excluded  from  the 
following  analysis  discussed  in  the  next  paragraph.  When  the 
sound  source  moves,  the  null  hypothesis  can  be  rejected  for 
all  the  twelve  points  tested  in  the  high-frequency  region. 
The  result  shows  that  for  high-frequency  signals,  the 
threshold  estimates  are  statistically  different  across  the 
five  test  sessions  when  the  sound  source  is  moved  during  the 
measurement . 

The  average  standard  deviation  across  the  two  frequency 
ranges  and  the  four  subjects  is  displayed  in  Fig.  3-1.  For 
the  1-kHz  signal,  the  average  standard  deviation  for  four 
subjects  was  essentially  the  same  for  both  conditions,  about 
0.7  dB . The  result  shows  that  the  threshold  variance  for  1- 
kHz  signal  is  about  0.7  dB  and  the  fitting  variance  is 
negligible.  The  movement  of  the  hard-wall  source  across  the 
five  sessions  does  not  change  the  variability  in  the  threshold 
estimation  for  the  1-kHz  signal. 

For  signals  above  7 kHz,  the  average  standard  deviation 
(excluding  the  three  points  that  did  achieve  significance 
without  movement,  see  Table  3-2)  is  0.55  dB  for  the  condition 
in  which  the  sound  source  does  not  move.  This  shows  that  the 
threshold  variance  for  high-frequency  signals  is  slightly 
smaller  than  that  of  the  1-kHz  signal.  This  may  be  due  to  the 
difference  in  the  slope  of  psychometric  function  at  difference 
frequencies.  This  point  will  be  discussed  in  greater  detail 


Standard  deviation 
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1 kHz  Above  7 kHz 


Figure  3-1.  The  average  standard  deviation 
subjects  at  1 kHz  and  frequencies  above  7 kHz. 


of 


the  four 
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in  a subsequent  experiment  (Exp.  3)  . The  average  standard 
deviation  for  high-frequency  signals  increases  to  2.9  dB  for 
the  condition  in  which  the  sound  source  is  intentionally 
moved.  This  result  implies  that  the  fitting  variance  for 
high-frequency  signals  is  much  larger  than  that  of  1-kHz 
signal.  The  movement  of  the  hard-wall  source  increases 
significantly  the  variability  in  the  estimates  of 
high-frequency  thresholds. 

The  conclusion  of  this  experiment  is  that  the  intra- 
subject threshold  variability  caused  by  the  movement  of  the 
hard-wall  source  is  much  greater  for  frequencies  above  7 kHz 
than  for  lower  frequencies. 

Experiment  2 


Introduction 

The  above  experiment  demonstrates  that  the  change  in 
sound  pressure  caused  by  the  hard-wall  movement  results  in  a 
variability  of  ± 2.9  dB  in  the  estimates  of  high-frequency 
thresholds.  Under  practical  test  circumstances,  headphones 
are  usually  used  as  the  sound  source.  The  movement  of 
headphones  during  replacement  should  cause  sound  pressure 
changes  in  the  ear  canal  and  would  have  a discernible  effect 
on  the  estimation  of  high-frequency  thresholds.  In  this 
experiment,  we  will  measure  the  change  of  sound  pressure  in 
the  ear  canal  due  to  headphone  placement  and  the  resulting 
fitting  variance  will  be  evaluated. 
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Method 
Subi ects 

Subjects  SI,  S2 , and  S3  from  the  experiment  1 
participated  in  the  experiment  2 . 

Stimuli 

It  was  found  out  that  the  headphones  employed  had  much 
worse  frequency  response  for  frequencies  over  10  kHz  than 
those  below  10  KHz.  A technique  called  "chirp"  signals  was 
employed  to  measure  the  sound  pressure  in  the  ear  canal 
(Neumann,  Uppenkamp,  & Kollmeier,  1993)  . Chirp  signals  rather 
than  the  "Golay  codes"  signals  were  used  because  chirp  signals 
provided  a higher  signal -to-noise  ratios  at  high  frequencies 
without  causing  distortions  to  the  headphone.  Chirp  signals 
are  sinusoidal  signals  with  an  instantaneous  frequency  that 
changes  monotonically  with  time.  The  spectral  energy  at  a 
certain  frequency  is  inversely  proportional  to  the  rate  of 
change  of  the  instantaneous  frequency.  Thus,  by  suitably 
selecting  the  time  course  of  the  instantaneous  frequency  of  a 
chirp  signal,  it  is  possible  to  manipulate  its  spectrum.  An 
important  characteristic  of  chirp  signals  is  that  they  have 
low  peak  factors  (defined  as  the  range  of  the  waveform 
amplitude  divided  by  its  root-mean-square  value) . Two 
different  chirp  signals  both  having  flat  frequency  spectrum 
were  generated:  one  for  the  frequency  range  from  0 to  10  kHz 
and  the  other  for  the  frequency  range  from  10  to  20  kHz. 
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Procedure 

The  sound  pressure  in  the  ear  canal  was  measured  by  a 
probe  tube  (8  french  in  diameter)  connected  to  a Knowles 
electret  microphone.  The  probe  tube  was  brought  into  the  ear 
canal  through  the  natural  valley  in  the  concha  just  under  the 
tragus.  A small  earmold  with  a tiny  hole  under  the  tragus  was 
used  to  hold  the  probe  tube  without  blocking  the  cavity  of 
concha.  The  outside  part  of  the  probe  tube  was  fixed  with  a 
surgical  tape  to  the  cheek  bone  during  the  measurements  and 
the  whole  probe  system  was  relative  immune  to  disturbance 
caused  by  headphone  placement . 

With  headphones  in  the  proper  position,  the  frequency 
responses  of  the  ear  canal  were  measured  using  the  two  chirp 
signals.  The  same  procedure  was  repeated  seven  times  while 
removing  and  replacing  the  headphones  after  each  measurement. 
There  was  no  intentional  displacement  of  the  headphones  during 
this  process.  Two  types  of  headphones  (superaural  Sennheiser 
HD  450  and  circumaural  Koss  HV  /l)  were  tested. 

Results  and  Discussion 

The  total  frequency  response  of  the  ear  canal  in  the 
frequency  range  from  0 to  20  khz  is  shown  in  Fig.  3-2.  It 
displays  the  average  frequency  response  and  the  standard 
deviation  across  seven  different  headphone  positions  for  one 
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Figure  3-2.  The  solid  curve  is  the  average  frequency  response 
in  the  ear  canal  of  subject  1 for  seven  headphone  positions. 
Senheiser  HD  450  headphones  were  the  transducers.  The  error 
bars  represent  plus  minus  one  standard  deviation  across  the 
seven  measurements  at  different  frequencies. 
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Table  3-3  Standard  deviation  of  sound  pressure  in  the  ear 
canal  across  seven  headphone  positions 


Subj  ects 

Headphone 

Mean  standard  deviat ion (dB) 

0-8  kHz 

8-16  kHz 

SI 

Sennheiser  HD  450 

1 . 93 

3 . 77 

Koss  HV  /I 

1.36 

3 . 03 

S2 

Sennheiser  HD  450 

1.31 

1 . 99 

Koss  HV  /I 

1.42 

2 . 12 

S3 

Sennheiser  HD  450 

1 . 93 

1 . 91 

Koss  HV  /I 

1 . 74 

3 . 17 

AVG 

1.62 

2 . 67 
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of  the  subjects.  In  audiometry,  the  frequency  range  from  125 
to  8 kHz,  is  considered  low  frequencies,  while  the  frequency 
range  from  8 to  16  kHz,  is  considered  as  high  frequencies. 
The  average  standard  deviations  for  these  two  frequency  ranges 
is  listed  in  Table  3-3.  They  indicate  the  expected  sound 
pressure  change  in  the  ear  canal  caused  to  headphone 
placement.  As  expected,  the  sound  pressure  change  in  the  ear 
canal  is  smaller  for  low-frequency  region  than  for  high- 
frequency  region.  The  average  standard  deviation  for  the 
three  subjects  and  two  types  of  headphones  is  shown  in  the 
bottom  row  of  the  table. 

The  average  standard  deviation  from  8 to  16  kHz  for  the 
three  subjects  and  two  types  of  headphones  is  approximately 
2.7  dB.  It  is  about  half  the  size  of  the  pressure  change  in 
the  ear  canal  resulting  from  the  hard-wall  source  movement 
(+5.0  dB,  as  estimated  in  Chapter  2) . Compared  with  the  hard- 
wall  source,  the  headphone  has  a smaller  acoustic  impedance. 
This  may  be  one  possible  reason  for  the  small  change  in  sound 
pressure  caused  by  headphone  placement.  This  result  implies, 
as  one  would  expect,  that  the  fitting  variance  for 
conventional  headphones  is  smaller  than  that  of  the  hard-wall 
source.  The  first  experiment  of  this  chapter  showed  that  the 
combination  of  threshold  variance  and  fitting  variance  for  the 
hard-wall  source  is  about  ± 2.9  dB  for  high  frequencies. 
Thus,  we  would  expect  that  the  variability  in  threshold 
estimates  caused  by  headphone  placement  does  not  exceed  this 
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value.  This  estimate  may  be  compared  with  values  obtained  in 
the  previous  literature.  Testing  fifty  subjects  using  a Koss 
(HV/X)  earphone,  Stelmachowicz  et  al . (1989)  showed  that  the 
standard  deviation  of  intra-subject  thresholds  was  about  2.9 
dB  from  8 to  16  kHz.  Using  the  high-frequency  audiometer 
developed  by  Stevens  et  al . (1987),  Green,  Kidd,  and  Stevens 
(1987) , averaging  about  forty  subjects,  estimated  the 
intra-subject  standard  deviation  to  be  about  1.7  dB . 
Stelmachowicz,  Beauchaine , Kalberer , Kelly,  and  Jesteadt 
(1988)  estimated  the  intra-subject  standard  deviation  to  be 
2.7  dB  using  the  same  audiometer.  Both  estimates  were  made 
over  the  frequency  range  from  8 to  14  kHz . 

In  conclusion,  the  sound  pressure  change  in  the  ear  canal 
resulting  from  different  positioning  of  headphones  is  smaller 
than  that  resulting  from  movement  in  the  position  of  a hard- 
wall  source.  The  influence  of  the  headphone  placement  on  the 
variability  of  threshold  estimates  should  be  smaller  than  + 3 
dB. 


Experiment  3 


Introduction 


There  were  two  different  reports  on  how  the  reliability 
of  threshold  estimates  depends  upon  test  frequency.  One 
report  showed  that  the  reliability  of  threshold  estimates 
deteriorated  as  a function  of  frequency  (for  example,  Green  et 
al.,  1987;  Stelmachowicz,  et  al . , 1988,  1989).  The  other 
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report  showed  that  the  reliability  of  high-frequency 
thresholds  was  as  good  as  that  of  low-frequency  thresholds 
(for  example,  Fletcher,  1965;  Dreisbach  & Frank,  1991; 
Dreschler  et  al . 1985)  . The  reliability  of  threshold 
estimates  mainly  depends  upon  fitting  variance  and  threshold 
variance  at  test  frequency.  The  last  experiment  showed  that, 
with  circumaural  or  supraaural  headphones  as  transducers,  the 
fitting  variance  for  high  frequencies  was  larger  than  low 
frequencies.  In  this  study,  we  investigated  one  possible 
factor  that  may  account  for  smaller  threshold  variance  for 
high  frequencies  than  for  low  frequencies. 

The  relationship  between  the  variability  of  the 
threshold  estimates  and  the  slope  of  the  psychometric  function 
can  be  expressed  as 

V=[p(l-p)]/(-^)2  Eq.  3-1 
dx 

where  V is  the  variance  of  the  stimulus  estimate,  dF/dx  is  the 
slope  or  derivative  of  the  psychometric  function  and  p(l-p)  is 
the  binomial  variance.  Taylor  (1971)  gives  the  formula  in  the 
first  equation  of  his  paper.  He  credits  the  formula  to 
Robbins-Monro  (1951)  and  Wetherill  (1966).  For  a 
psychophysical  procedure,  we  track  a specific  percent -correct 
point  (for  example,  p = 0.79  for  3 -down  1-up  2 AFC) . Thus  the 
binomial  variance,  p(l-p)  is  the  same  at  the  threshold  point 
for  different  stimuli.  Therefore,  the  variability  of  the 
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threshold  estimates  is  inversely  proportional  to  the  slope  of 
the  psychometric  function. 

Watson,  Franks,  and  Hood  (1972)  found  that  the  slope  of 
the  psychometric  function  for  pure-tone  signals  increases  from 
125  to  4000  Hz.  If  this  holds  for  frequencies  above  8 kHz, 
then  the  steeper  psychometric  functions  at  high  frequencies 
would  reduce  the  variability  in  the  estimates  of  the  absolute 
thresholds.  In  this  experiment,  we  measured  the  slope  of  the 
psychometric  function  at  1 and  15  kHz. 

Method 

Subi ects 

Four  subjects  with  normal  hearing  sensitivity  (250  to 
8000  Hz)  participated  in  the  experiment.  All  subjects  were 
students  at  the  University  of  Florida.  The  average  age  of  the 
subjects  was  22  years.  An  earmold  was  made  for  each  subject 
and  the  subjects  were  asked  to  have  their  external  ears 
cleaned  before  the  experiment.  Prior  to  data  collection,  each 
subject  completed  several  hours  of  practice  on  the  task. 
Subjects  were  paid  an  hourly  wage  and  received  a 20  percent 
bonus  upon  completion  of  the  study. 

Stimuli  and  procedure 

An  adaptive  procedure  was  used  to  measure  the 

psychometric  function  at  1 and  15  kHz  (Dai,  1994) . The  level 
of  the  pure-tone  signal  was  changed  adaptively  following  the 
three-down,  one-up  method  (Levitt,  1971)  in  a two-interval. 
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forced-choice  procedure  (2AFC) . The  step  size  was  always 
fixed  at  1 dB  throughout  each  block  of  trials  and  the  signal 
had  a duration  of  200  ms.  This  procedure  has  a major 
advantage  over  the  conventional,  constant -stimulus  procedure. 
The  constant -stimulus  procedure  requires  that,  prior  to 
measuring  the  psychometric  functions,  the  experimenter  has  to 
determine  signal  levels  that  will  cover  a range  of  response 
between  about  55%  to  95%  correct.  In  contrast,  owing  to  its 
adaptive  nature,  the  present  procedure  automatically  adjusts 
the  signal  levels  to  be  within  the  desired  range  of  percent- 
correct  values. 

Leek,  Hanna,  and  Marshall  (1992)  examined  the  reliability 
of  psychometric  functions  measured  using  an  up-down  procedure. 
Based  on  simulations,  they  concluded  that  the  adaptive 
procedure  gave  reasonably  accurate  estimate  of  the  true 
psychometric  functions  when  at  least  200  trials  were  included 
in  obtaining  each  function.  In  the  present  study,  each 
psychometric  function  was  based  on  a block  of  600  trials.  For 
each  signal  frequency,  8 blocks  of  data  were  obtained. 

Results  and  Discussion 

In  fitting  the  psychometric  functions,  we  only  used  data 
that  met  the  following  two  criteria.  First,  only  signal 
levels  that  received  more  than  20  trials  were  kept.  Second, 
signal  levels  that  yielded  percent -correct  scores  between  one 
standard  deviation  ( [p  (1-p) /n] 1/2  ) below  100%  and  one 
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Table  3-4.  An  example  of  psychometric  function  at  1 kHz 


Signal  level 

# of  Trials 

# correct 

% correct 

0 

9 

9 

100 

-1 

12 

12 

100 

-2 

22 

21 

95 

-3 

59 

55 

93 

-4 

114 

98 

86 

-5 

155 

130 

84 

-6 

131 

98 

75 

-7 

73 

49 

67 

-8 

22 

13 

59 

-9 

3 

1 

33 

Total 

600 

486 
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standard  deviation  above  50%  were  selected  to  fit  the 
psychometric  functions.  Using  these  criteria,  about  8%  of  the 
data  were  excluded  from  the  analysis  of  the  psychometric 
functions.  Table  3-4  shows  an  example  of  a single 
psychometric  function  which  was  measured  at  1 kHz  for  subject 
S3  . In  fitting  the  psychometric  function,  we  only  used  signal 
levels  between  -8  and  -2  dB . 

Egan,  Lindner,  and  McFadden  (1966)  have  shown  that  the 
relation  between  detection  scores  in  the  2AFC  procedure  and 
the  ratio  of  signal  energy  to  noise  power  per  unit  bandwidth 
(E/NO)  is  reasonably  well  described  by 

K Eq.  3-2 
No 

Since  both  signal  duration  and  noise  spectrum  level  could  be 
treated  as  constants  in  the  present  experiment,  Eq.  3-2  could 
be  simplified  to  a logarithmic  version 

lOlog  (d')=KS+B  Eq.  3-3 

where  S is  the  sound  pressure  level  of  the  signal  (in 
decibels),  and  K and  B are  fitting  constants.  In  the  fitting 
procedure,  the  detectability  index  d'  was  transformed  to 
percent-correct  scores  (Elliot,  1964;  Green  & Dai,  1991)  . 

The  error  to  be  minimized  is 
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PEi-PMi 
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Eq.  3-4 


where  PEi  and  PMi  are  expected  and  measured  percent -correct 
scores,  cq  is  the  standard  deviation  of  PMi7  and  N is  the 
number  of  points  to  be  fitted.  K and  B for  the  best  fit  were 
obtained . 

The  slopes  of  the  fitted  psychometric  functions,  K,  for 
three  subjects  are  shown  in  Table  3-5.  Comparing  the  average 
slopes  of  the  fitted  psychometric  functions  for  different 
frequencies,  the  slope  at  15  kHz  is  steeper  than  that  at  1 kHz 
for  all  three  subjects.  This  is  consistent  with  the  results 
reported  by  Watson  et  al . (1972)  . They  showed  that  the  fitted 
psychometric  functions  were  steeper  at  high  frequencies  than 
at  low  frequencies  from  125  to  4000  Hz.  Combining  data  for 
three  subjects,  a one-way  ANOVA  also  supports  the  conclusion 
that  the  slope  of  the  psychometric  function  is  greater  at  15 
kHz  than  at  1 kHz  (p  < 0.05) . The  average  increase  in  slope 
from  1 to  15  kHz  is  only  about  25  percent.  According  to  Eq. 
3-1,  the  increase  in  the  slope  of  the  psychometric  function 
may  result  in  25%  smaller  standard  deviation  in  the  estimates 
of  absolute  threshold.  This  agrees  with  the  result  on 
threshold  variance  in  Experiment  1. 

The  mean  slope  of  the  fitted  psychometric  function  at 
1000  Hz  is  0.97  in  this  study.  This  corresponds  to  a range  of 
6-10  dB  from  near-chance  (55%)  to  near-perfect  (95%) 
performance.  It  agrees  with  the  results  obtained  in  several 


Table  3-5.  Slopes  of  fitted  psychometric 
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previous  studies  using  various  psychophysical  procedures  (Egan 
et  al . , 1966;  Green  & Swets,  1966;  Watson,  Franks,  & Hood, 
1967;  Barr-Hamilton,  Bryan,  & Tempest,  1968)  . Watson  et  al . 
(1972),  however,  obtained  much  steeper  psychometric  functions 
in  their  study.  Their  calculated  psychometric  function  had  a 
range  of  3-6  dB . They  indicated  that  if  the  psychometric 
function  was  calculated  based  on  data  pooled  over  all  test 
sessions,  the  obtained  slope  was  approximately  1.  The  steeper 
slope  was  obtained  when  the  psychometric  function  was 
calculated  based  on  data  in  isolated  sessions,  in  which  on 
headphone  replacement  occurred.  In  our  study,  each  value  of 
the  slope  was  calculated  from  the  data  in  each  test  session  in 
which  the  experimental  condition  remained  the  same.  The 
discrepancy  between  the  results  obtained  by  Watson  et  al . 
(1972)  and  other  studies  is  not  clear. 

In  conclusion,  for  all  the  three  subjects  tested,  the 
psychometric  function  at  15  kHz  is  steeper  than  that  at  1 kHz. 
This  may  contribute  to  the  low  variability  in  the  estimates  of 
high-frequency  thresholds. 


Summary 

In  this  chapter,  we  investigated  the  reliability  of 
intra-subject  threshold  caused  by  transducer  placement.  The 
first  transducer  being  investigated,  the  hard-wall  source,  had 
very  high  acoustic  impedance  and  would  create  very  strong 
standing  waves  in  the  ear  canal . The  results  from  the 
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previous  chapter  showed  that  the  change  in  sound  pressure  at 
the  ear  canal  was  approximately  ± 5 dB  if  the  sound  source 
moved  by  one-quarter  wavelength.  The  variation  in  behavior 
threshold  was  0.7  dB  for  1-kHz  signal  and  2.9  dB  for  signals 
above  7 kHz  if  the  sound  source  moved  by  as  much  as  ± 0.5  cm. 
The  smaller  variation  in  low-frequency  threshold  is  due  to  its 
larger  wavelength.  It  was  found  out  that  the  fitting  variance 
for  high-frequency  signals  were  much  larger  than  that  of  the 
low-frequency  signals. 

Next,  two  types  headphones  frequently  used  in  clinics 
were  investigated.  The  sound  pressure  change  in  the  ear  canal 
caused  by  headphone  placement  was  found  to  be  smaller  than 
that  caused  by  hard-wall  source  movement.  This  is  probably 
due  to  the  smaller  acoustic  impedance  for  the  headphones. 
Thus  for  high-frequency  signals,  the  resulting  fitting 
variance  measured  with  the  headphones  was  expected  to  be 
smaller  than  ± 2.9  dB  as  measured  with  high-wall  source. 

Finally,  the  slopes  of  the  psychometric  functions  for  1- 
kHz  and  15 -kHz  signals  were  measured.  It  was  found  that  the 
slope  of  the  psychometric  function  for  15 -kHz  signal  is  25% 
larger  than  that  of  the  1-kHz  signal.  This  may  result  in  25% 
smaller  standard  deviation  in  the  estimates  of  absolute 
threshold  for  high-frequency  than  low-frequency  signals.  The 
larger  fitting  variance  and  smaller  threshold  variance  for 
high-frequency  signals  may  be  responsible  for  the  diverse 
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reports  on  the  reliability  of  high-frequency  threshold 
estimates . 


Chapter  4 

ESTIMATES  OF  FREQUENCY  RESOLUTION  AT  HIGH  FREQUENCIES 

Introduction 

The  analytic  property  of  the  auditory  system,  the  ability 
to  hear  out  individual  frequency  components  of  a complex  sound 
was  first  described  in  Ohm's  (1843)  acoustic  law.  Later 
studies  by  Plomp  (1964)  and  Duifhuis  (1970,  1971)  have  refined 
the  earlier  estimates.  Bekesy' s investigations  (1960) 
provided  direct  observation  of  the  filtering  action  of  the 
cochlea  and  it  is  generally  conceded  that  the  ear  can  be 
treated  (modeled)  as  a bank  of  overlapping  bandpass  filters. 

Psychoacoustic  experiments  frequently  used  masking  as  a 
way  to  determine  the  degree  and  precision  of  the  filtering 
action.  In  a masking  experiment,  one  sound  is  made  more 
difficult  to  detect  by  the  presence  of  another  (masking) 
sound.  Wegel  and  Lane  (1924)  studied  how  one  sinusoid  masked 
another  and  presented  dramatic  evidence  that  a listener  could 
attend  to  a single  frequency  and  ignore  a wide  range  of 
masking  energy  unless  its  frequency  was  very  close  to  the 
signal  being  attended. 

Fletcher's  (1940)  classical  experiment  on  the  masking 
produced  by  noise  laid  the  foundation  for  the  concept  of  what 
we  now  call  the  critical  band.  In  Fletcher's  experiment,  the 
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detectability  of  a sinusoidal  signal  centered  within  a band  of 
(masking)  noise  was  estimated  as  a function  of  the  bandwidth 
of  the  noise.  Starting  with  a very  narrow  band  of  noise,  the 
threshold  of  the  signal  increased  as  the  noise  bandwidth 
increased.  Then  above  some  critical  bandwidth  point,  which 
depended  on  the  signal  frequency,  further  increases  in  the 
noise  bandwidth  did  not  result  in  significant  changes  in 
signal  threshold.  Fletcher  argued  that  some  kind  of  "internal 
filter"  was  centered  at  the  signal  frequency  and  that  the 
detection  of  the  signal  was  determined  by  the  ratio  of  signal 
and  noise  power  passed  by  this  filter.  When  the  noise 
bandwidth  increased  initially,  the  noise  power  passed  through 
the  internal  filter,  often  termed  an  auditory  filter, 
increased.  When  the  noise  bandwidth  was  greater  than  the 
auditory  filter  bandwidth,  further  widening  in  the  noise 
bandwidth  did  not  increase  the  noise  power  passing  through  the 
auditory  filter.  Fletcher  called  the  bandwidth  at  which  the 
signal  threshold  ceased  to  increase  the  "critical  bandwidth" . 

Since  Fletcher's  experiment,  many  measurements  of  the 
auditory  critical  bands  have  been  undertaken  (Hamilton,  1957; 
Swets,  Green,  & Tanner,  1962;  Houtgast,  1974;  Patterson,  1976; 
Patterson  & Henning  1977) . Practically  no  measurements  have 
been  obtained  at  frequencies  above  10  kHz.  One  of  the 
problems  is  that  it  is  difficult  to  generate  these  high 
frequency  stimuli  with  controlled  spectra  using  traditional 
headphones.  Another  problem  is  the  variability  in  the  high- 
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frequency  threshold  estimates  introduced  by  the  fitting 
variance  as  discussed  in  the  previous  chapter.  These  two 
problems  can  be  avoided  by  using  the  sound  delivery  system 
described  in  Chapter  2.  In  this  chapter,  two  different 
methods  were  used  to  estimate  auditory  frequency  resolution  at 
frequencies  above  10  kHz.  One  method  is  the  modification  of 
Fletcher's  original  masking  experiment.  The  second  approach 
is  based  on  the  monaural  phase  discrimination  of  a complex 
tone . 

Experiment  1 

Introduction 

In  this  study,  the  auditory  filter  shapes  were  estimated 
using  the  well-established  masking  method  in  which  a notched- 
noise  masker  was  employed  (Patterson,  1976;  Patterson  & Nimmo- 
Smith,  1980;  Patterson,  Nimmo-Smith,  Weber,  & Milroy,  1982; 
Glasberg  & Moore,  1986;  Glasberg  & Moore,  1990,  Shailer, 
Moore,  Glasberg,  Watson,  & Harris,  1990)  . Moore  and  Glasberg 
(1990)  summarized  the  data  on  the  estimates  of  the  equivalent 
rectangular  bandwidth  (ERB)  of  the  auditory  filter  for  center 
frequencies  under  10  kHz.  There  was  no  measurement,  however, 
for  center  frequencies  over  10  kHz . 

The  knowledge  of  the  auditory  filter  shapes  above  10  kHz 
is  important  in  the  understanding  of  the  frequency  selectivity 
of  the  auditory  system  in  several  aspects.  A number  of 
experiments  have  showed  that  high  frequency  hearing  is  first 
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and  most  damaged  by  noise  or  ototoxic  drugs  (Fletcher  et  al . , 
1967;  Fausti,  et . al . 1981;  Jacobson,  et  al . , 1965).  With  the 
knowledge  of  the  auditory  filter  shape,  the  masking  caused  by 
noise  energy  in  the  high-frequency  region  can  be 
quantitatively  estimated  through  the  calculation  of  masking 
and  excitation  patterns.  The  influence  of  noise  and  ototoxic 
drugs  on  the  frequency  selectivity  of  the  patients  can  also  be 
examined  in  the  high-frequency  regions.  Evidence  showed  that 
energy  at  high  frequencies  and  spectral  shape  cues  were  very 
important  to  the  accurate  vertical  localization  (Roffler  & 
Butler,  1968;  Gardner  & Gardner  1973;  Butler  & Helwig,  1983; 
Middlebrooks  & Green,  1990)  . The  knowledge  of  the  auditory 
filter  shapes  at  high-frequencies  may  help  us  analyze  the 
spectral  shape  cues  in  the  localization  of  sound  source 
elevation . 

Method 

Subi ects 

All  three  subjects  were  students  at  the  University  of 
Florida.  Two  of  them  were  recruited  from  an  advertisement  in 
the  campus  newspaper  and  were  paid  for  their  participation. 
The  other  one  was  the  author  of  the  thesis.  The  average  age 
of  the  subjects  was  21.  Each  subject's  audiogram  when 
measured  with  a Bekesy  audiometer  was  within  10  dB  of  normal 
hearing  (ISO  standard)  in  the  range  of  125  to  8000  Hz.  All 
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three  subjects  had  previous  experience  in  psychoacoustic 
experiments . 

Stimuli 

The  sinusoidal  signal  and  the  masking  noise  were 
generated  digitally  using  an  array  processor  (TDT  QAP1) 
controlled  by  a IBM  PC  compatible  computer.  The  sinusoidal 
signal  and  the  masking  noise  were  delivered  via  two  separate 
channels  of  a digital -to-analog  converter  (TDT  QDA2 ) and  their 
levels  were  controlled  by  two  programmable  attenuators.  The 
signal  and  masker  were  combined  in  an  adder  before  being 
passed  to  the  transducer  (Radio  shack  40-1310B)  through  a 
power  amplifier  (Crown  D60)  . The  signal  generated  by  the 
transducer  was  delivered  to  the  ear  canal  through  the  sound 
delivery  system  described  in  Chapter  2. 

In  order  to  generate  a noise  masker  with  a flat  spectrum 
in  the  desired  frequency  range  to  the  ear  canal,  the  input 
electrical  signal  to  the  transducer  was  inverse  filtered  by 
the  frequency  response  of  the  sound  delivery  system.  The 
pressure  delivered  to  the  ear  canal  was  measured  with  an 
electret  microphone  (Knowles  Electronics,  EA-1954)  mounted  at 
the  end  of  the  long  lossy  tube.  The  spectrum  of  the  stimuli 
was  monitored  with  a spectrum  analyzer  (Data  Precision,  Data 
6000) . The  recorded  spectrum  of  a Gaussian  noise  with  a flat 
spectrum  from  1 to  19  kHz  is  shown  in  Fig.  4-1.  It  shows  that 
we  can  control  the  spectrum  of  the  stimuli  delivered  to  the 
ear  canal  with  good  accuracy.  The  Knowles  electret  microphone 
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was  calibrated  against  a standard  Briiel  & Kjaer  microphone 
(model  number  4133)  . 

The  notched-noise  masker  was  generated  digitally  in  the 
present  study,  which  is  different  from  the  modulation 
technique  described  by  Patterson  (1976)  . Broad  band  Gaussian 
noise  was  generated  in  the  frequency  domain  by  filling  a 8192 
point  buffer  with  numbers  having  Gaussian  distribution  and 
zero  mean.  The  notch  was  generated  by  setting  the  number 
outside  the  desired  passband  to  zeros.  The  notch  noise 
waveform  was  obtained  by  taking  the  inverse  Fourier 
transformation  of  the  defined  spectrum. 

Previous  reports  showed  that  it  was  more  difficult  to 
detect  the  signal  from  the  masker  when  the  signal  and  the 
masker  were  gated  on  and  off  synchronously  (Scholl,  1962; 
Carlyon,  1989;  Shailer  et  al . , 1990).  In  this  study,  the 
signal  was  switched  on  150  ms  after  the  onset  of  the  noise  and 
switched  off  50  ms  before  the  offset  of  the  noise.  The 
masking  noise  and  the  signal  are  400  ms  and  200  ms  in  duration 
respectively.  Both  the  masker  and  the  signal  had  a 10  ms 
rise/fall  time.  The  signal  frequencies  tested  were  12,  and  14 
kHz  for  the  high-frequency  region.  In  order  to  make  a 
comparison  with  previous  experiments  measured  with  headphone 
systems,  we  also  tested  4 kHz. 

The  notch  width  of  the  noise  masker  is  specified  as  the 
deviation  of  the  nearest  edges  of  the  lower  and  upper  noise 
bands  from  center  frequency  fc,  divided  by  the  center 
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Figiire  4 1.  Recorded  spectrum  of  a Gaussian  noise  at  the  input 
to  the  ear  canal . The  Gaussian  noise  had  a flat  spectrum  from 
1 to  19  kHz. 
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frequency.  This  will  be  denoted  as  Af/fc.  The  experiment  had 
three  basic  sets  of  conditions.  For  the  first  set  of 
condition,  the  notch  was  placed  symmetrically  around  fc.  The 
values  of  notch  width,  expressed  as  Af/fc,  were  0.0,  0.1,  0.2, 
0.3,  0.4,  0.5,  0.6,  and  0.7.  For  the  second  set  of  condition, 
the  near  edge  of  the  upper  band  was  0.2  further  away  from  fc 
than  the  near  edge  of  the  lower  band.  The  values  of  Af/fc  for 
the  lower/upper  bands  were:  0.1/0.3,  0.2/0. 4,  0.3/0. 5, 
0.4/0. 6,  and  0.5/0. 7.  The  third  set  of  conditions  was  a 
mirror  image  of  the  second  set  of  conditions,  with  the  lower 
band  being  further  away  from  fc  than  the  upper  band.  The 
values  of  Af/fc  were:  0. 3/0.1,  0.4/0. 2,  0.5/0. 3,  0.6/0. 4,  and 
0.7/0. 5.  Each  noise  band  had  a width  of  0.4  fc.  At  the 
widest  notch  widths,  the  width  of  the  lower  and  upper 
frequency  bands  were  adjusted  to  avoid  " wrapping  around"  zero 
frequency  and  Nyquist  frequency  (20,400  Hz). 

Procedure 

A two-interval  forced-choice  adaptive  procedure  (Levitt, 
1971)  was  used  to  determine  signal  thresholds  in  blocks  of  60 
trials.  A trial  consists  of  two  400  ms  observation  intervals 
with  a 500  ms  inter-stimulus  silent  interval  between  them. 
Each  block  began  with  the  signal  about  10  dB  above  the  masked 
threshold.  After  three  consecutive  correct  responses,  the 
signal  level  was  decreased  and  for  each  incorrect  response, 
the  signal  level  was  increased.  The  step  size  was  5 dB  at  the 
beginning  of  the  block,  but  after  two  reversals,  it  was 
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decreased  to  2 dB . This  procedure  estimates  the  79.4  % point 
on  the  psychometric  function.  Each  threshold  point  was  the 
average  of  5 blocks . 

The  absolute  thresholds  of  the  subjects  were  measured  at 
8,  10,  12,  14,  16  kHz  using  the  same  2AFC  procedure. 

To  insure  that  the  position  of  the  source  did  not  change 
during  the  course  of  the  experiment,  the  impulse  response  of 
the  ear  canal  was  recorded  at  the  beginning  and  end  of  each 
notch  width  condition.  For  all  the  data  reported,  the  rms 
difference  between  the  two  impulse  responses  was  less  than  1.5 
dB  in  the  frequency  range  from  5000  to  16000  Hz.  The  earmold 
was  removed  after  each  condition,  and  the  subject  took  a ten- 
minute  break.  At  the  beginning  of  a second  condition,  the 
positions  of  the  earmold  and  the  plastic  tube  were  adjusted  so 
that  the  rms  of  the  impulse  response  was  within  1.5  dB  of  that 
recorded  at  the  end  of  the  previous  condition. 

Calculation  of  the  Shapes  of  the  Auditory  Filters 

The  procedure  of  deriving  the  auditory  filter  shapes  is 
very  similar  to  that  described  by  Glasberg  and  Moore  (1990) . 
A brief  description  is  given  here. 

According  to  the  power- spectrum  model  of  masking,  the 
threshold  Ps  for  a signal  in  noise  is 


Eq.  4-1 
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where  K is  a constant,  f is  frequency,  N(f)  represents  the 
spectrum  of  noise  masker  at  the  input  to  the  auditory  filter 
and  W(f)  represents  the  intensity  weighting  function  of  the 
auditory  filter. 

In  this  study,  the  function  N(f)  is  evaluated  in  two 
different  methods.  For  the  first  method,  the  input  noise 
spectrum  to  the  auditory  filter  is  considered  to  be  flat 
within  its  passbands . The  signal  threshold  can  be  predicted 
by  integrating  the  intensity  weighting  function  ranged  within 
the  noise  passband.  This  can  be  expressed  as 

fl/fo  Afu/f0  + 0.4 

Ps=KN0f0  j*  W(g)dg+kN0f0  J W{g)  dg  Eq.  4-2 

-0.4-Afj/fo  a fu/f0 

where  Aflf  Afu,  are  the  lower  and  upper  notch  width,  g is  a 
normalized  frequency  variable  (g  = |f-fc|/f0)  and  N0  is  the 
spectrum  level  of  the  noise  masker. 

For  the  second  method,  the  input  noise  spectrum  to  the 
auditory  filter  is  corrected  for  the  filtering  function  of  the 
external  and  middle  ear.  Because  the  critical-band  filtering 
occurs  after  the  middle  ear,  a frequency-dependent 
attenuation  function  caused  by  the  transfer  function  of  the 
external  and  middle  ear  is  applied  to  the  flat  noise  masker. 

According  to  Dallos  (1973),  the  transducers  in  the 
cochlear  are  equally  sensitive  to  stimuli  across  all  audible 
frequencies . Thus  the  absolute  threshold  curve  across 
frequency  reflects  the  frequency-dependent  function  resulting 
from  the  transfer  function  of  the  external  and  middle  ear. 
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Following  Glasberg  and  Moore  (1990) , for  frequencies  below  8 
kHz,  we  use  the  minimum  audible  pressures  (MAP)  estimated  by 
Killion  (1978)  as  the  absolute  threshold  curve.  For 
frequencies  above  8 kHz,  the  measured  average  absolute 
threshold  for  the  three  subjects  are  used.  The  absolute 
threshold  curve  is  extrapolated  to  20  kHz  using  a polynomial 
function.  The  measured  absolute  thresholds  for  the  four 
subjects  are  shown  in  Table  4-1. 

The  filter  weighting  function  W(g)  is  assumed  to  have  a 
general  form  of  a round  exponential  function  (roex) . Two 
different  roex  functions  were  used  to  fit  the  data.  One  is 
known  as  roex(p,r)  model.  In  this  model,  each  side  of  the 
auditory  filter  is  assumed  to  have  the  form 

W(g)  = (1-r)  (1+pg)  exp(-pg) +r  Eq.  4-3 

where  p is  a parameter  determining  the  slope  of  the  filter 
skirts,  and  r is  a parameter  that  flattens  the  tails  of  the 
filter . 

The  other  model  is  known  as  roex(p,w,t)  model.  In  this 
model,  each  side  of  the  auditory  filter  is  assumed  to  have  the 
form 

W{g)  = (1-wj  (l+pgr)  exp  {-pg)  + w(l+tg)  exp  (-tg)  Eq.  4-4 

where  p is  a parameter  determining  the  slope  of  the  filter 
skirts,  and  t is  a parameter  determining  the  shallower  tail  of 
the  filter,  and  w determines  where  the  tail  takes  over  from 
the  passband. 
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Table  4-1:  Absolute  thresholds  for  the  individual  subjects  at 
frequencies  of  8,  10,  12,  14  and  16  kHz,  in  dB  SPL.  The  means 
of  the  three  subjects  were  used  to  define  the  MAP  correction 


Subj  ect 

8 

10 

Frequency, 

12 

kHz 

14 

16 

SI 

13 . 0 

20 . 0 

16 . 0 

14 . 0 

26 . 0 

S2 

19.3 

27.8 

16 . 6 

14 . 3 

36 . 7 

S3 

16 . 7 

10 . 9 

13.8 

12 . 9 

42 . 6 

Mean 

16.3 

19 . 6 

15.5 

13 . 7 

35 . 1 
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We  found  out  the  roex(p,w,t)  model  fit  the  data  better 
than  the  roex(p,r)  model.  The  subsequent  analysis  will 
concentrate  on  the  results  obtained  with  the  roex(p,w,t) 
model.  For  the  roex(p,w,t)  model,  without  any  restriction,  it 
has  six  free  parameters,  three  for  each  side  of  the  filter. 
According  to  Paterson  and  Nimmo-Smith  (1980)  and  Glasberg, 
Moore,  Patterson,  and  Nimmo-Smith  (1984),  the  model  can  be 
simplified  by  letting  wu  = wx  = w,  and  tu/tx  = pu/p1#  where  the 
subscript  u and  1 represent  the  upper  and  lower  side  of  the 
filer.  With  the  simplification,  the  number  of  free  parameters 
was  reduced  to  four. 

The  fitting  procedure  works  in  the  following  iterative 
process.  Starting  with  a set  of  values  of  plf  pu,  and  w and 
tx.  For  each  notch  width,  the  center  of  the  auditory  filter 
is  shifted  to  the  position  that  gives  the  highest  signal-to- 
noise  ratio.  Equation  4-1  is  then  used  to  predict  the 
threshold  for  each  notch  width  condition.  The  mean- squared 
deviation  between  the  measured  thresholds  and  the  predicted 
threshold  is  minimized  by  changing  the  values  of  Plf  Pu,  and 
w and  t1 . 

Results  and  Discussion 

The  individual  and  mean  threshold  data  at  different  notch 
width  conditions  are  shown  by  three  symbols  in  Figs.  4-2  to  4- 
6.  The  open  circles  denote  conditions  that  the  noise  bands 
were  symmetrically  placed  about  the  signal  frequency.  The 
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right  arrows  denote  conditions  that  the  upper  noise  bands  were 
further  away  from  the  center  frequency  than  the  low  noise 
bands.  The  left  arrows  denote  the  opposite  conditions.  For 
a center  frequency  of  4 kHz,  the  noise  spectrum  level  was  40 
dB  SPL.  For  center  frequencies  of  12  and  14  kHz,  noise 
spectrum  levels  of  4 0 and  50  dB  were  tested.  The  solid  curves 
in  each  panel  give  the  thresholds  fitted  to  the  measured  data 
using  the  roex(p,w,t)  model  without  MAP  correction.  Figures 
4-7  to  4-9  show  the  filter  shapes  derived  from  the  threshold 
curves  at  different  center  frequencies  and  spectrum  levels. 
Table  4-2  list  the  parameters  for  the  auditory  filters  fitted 
to  the  measured  data  both  with  and  without  the  MAP  correction. 
It  also  gives  the  rms  deviation  of  the  data  from  the  fitted 
values  for  each  subject  and  condition  and  the  value  of  K in 
Eq.  4-1. 

The  effect  of  including  the  MAP  correction 

Contrary  to  the  results  reported  by  Shailer  et  al . 
(1990)  , we  found  that  the  fit  to  data  was  not  better  with  the 
MAP  correction.  The  root -mean- square  deviation  of  the  data 
from  the  fitted  values  is  smaller  with  the  MAP  correction  than 
without  the  correction  for  only  six  of  the  twenty  fitted 
filters.  Thus  we  will  discussed  the  filter  shapes  derived 
without  the  MAP  correction  in  the  rest  of  the  paper.  For 
completeness,  the  parameters  of  the  filters  derived  with  the 
MAP  correction  are  also  given  in  Table  4-2. 


f0  = 4 kHz,  N0  = 40  dB 
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DEVIATION  OF  NEARER  NOISE  EDGE.  Af/f 


Figure  4-2.  Individual  and  mean  threshold  data  for  a center 
frequency  of  4 kHz  and  a spectrum  level  of  40  dB . Asterisks 
denote  conditions  where  the  noisebands  were  symmetrically 
placed  about  the  signal  frequency;  right -pointing  arrows 
indicate  that  the  notch  was  centered  above  the  signal 
frequency,  and  left-pointing  arrows  indicate  that  the  notch 
was  centered  below  the  signal  frequency.  The  abscissa  refers 
to  the  distance  between  the  nearer  edge  of  the  two  noise  bands 
to  the  signal  frequency.  The  solid  lines  are  the  thresholds 
fitted  the  data  by  the  roex(p,w,t)  model. 


f0  = 12  kHz,  N0  = 40  dB 
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DEVIATION  OF  NEARER  NOISE  EDGE,  Af/f 


Figure  4-3.  As  in  Fig.  4-2  but  for  a center  frequency  of  12 
kHz  and  a spectrum  level  of  40  dB . 


f0  = 12  kHz.  N0  = 50  dB 
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DEVIATION  OF  NEARER  NOISE  EDGE.  Af/f 


Figure  4-4.  As  in  Fig.  4-2  but  for  a center  frequency  of  12 
kHz  and  a spectrum  level  of  50  dB . 


f0  = 14  kHz.  N0  = 40  dB 
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DEVIATION  OF  NEARER  NOISE  EDGE,  Af/f 


Figure  4-5.  As  in  Fig.  4-2  but  for  a center  frequency  of  14 
kHz  and  a spectrum  level  of  40  dB . 


f0  = 14  kHz,  N0  = 50  dB 
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DEVIATION  OF  NEARER  NOISE  EDGE,  Af/f 


Figure  4-6.  As  in  Fig.  4-2  but  for  a center  frequency  of  14 
kHz  and  a spectrum  level  of  50  dB . 
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f0  = 4 kHz.  N0  = 40  dB 


DEVIATION  FROM  CENTER  FREQUENCY,  Af/f 


Figure  4-7.  Filter  shapes  derived  from  the  individual  and  the 
mean  data  for  a center  frequency  of  4 kHz  and  a spectrum  level 
of  40  dB . 
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DEVIATION  FROM  CENTER  FREQUENCY.  Af/f 


Figure  4-8.  Filter  shapes  derived  from  the  individual  and  the 
mean  data  for  a center  frequency  of  12  kHz  and  spectrum  levels 
of  40  dB  (solid  line)  and  50  dB  (dotted  line) . 
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f„  = 14  kHz 


DEVIATION  FROM  CENTER  FREQUENCY.  Af/f 


Figure  4-9.  As  in  Fig.  4-8  but  for  a center  frequency  of  14 
kHz . 


Table  4-2.  Summary  of  the  parameters  of  the  auditory  filters  fitted  to  the  data.  Table 
also  shows  the  values  of  K,  as  in  Eq.  4-1,  and  the  root -mean- square  (rms)  deviation  of 
the  data  from  the  fitted  values.  The  upper  number  of  each  pair  was  obtained  without  any 
correction;  the  number  in  parentheses  was  obtained  using  the  MAP  correction.  For  each 
frequency  and  level,  the  first  mean  is  the  mean  value  of  the  parameters  across  subjects. 
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The  shapes  of  the  derived  filters 

At  the  low  center  frequency  of  4 kHz,  for  all  three 
subjects,  the  slope  of  the  filter  is  shallower  at  low- 
frequency  side  than  at  high-frequency  side.  This  is 
consistent  with  results  reported  in  previous  studies  (for 
example,  Glasberg  et  al . , 1984) . 

At  center  frequencies  of  12  and  14  kHz,  the  filters 
become  broader  for  50  dB  level  than  for  40  dB . This  trend  is 
consistent  with  data  at  low  center  frequencies  showing  that 
the  auditory  filter  becomes  broader  with  increasing  masking 
level  (Weber,  1977) . On  the  other  hand,  although  earlier  work 
at  lower  center  frequencies  suggested  that  the  high-frequency 
side  of  the  filter  changes  only  slightly  with  level,  or  even 
increases  in  sharpness  with  increasing  level  (Lutfi  & 
Patterson,  1984;  Moore  & Glasberg,  1987),  the  present  results 
show  a broadening  on  the  high-frequency  side  at  higher 
spectrum  level . 

It  is  interesting  to  notice  that  for  subject  1,  he  has 
the  broadest  auditory  filters  for  all  the  center  frequencies 
and  spectrum  levels.  The  bandwidths  of  auditory  filters  for 
subject  2 and  3 are  similar. 

The  goodness  of  fitting  with  the  roex(p.w.t)  model 

Using  the  roex(p,w,t)  model,  it  was  found  that  the  fit  to 
the  data  was  much  better  for  the  lower  spectrum  level  than  the 
higher  spectrum  level.  At  40  dB  SPL,  the  mean  rms  values  for 
4,  12,  and  14  kHz  are  2.1,  1.6  and  2.8  dB,  respectively.  At 
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50  dB  SPL,  the  mean  rms  values  for  12  and  14  kHz  are  2.7  and 
4.3  dB,  respectively.  For  several  fitted  filters,  especially 
at  50  dB  SPL,  the  fitted  parameters  tx  and  tu  are  zero.  This 
means  that  the  tail  part  of  the  filter  does  not  contribute  to 
the  overall  filter  shape,  and  the  effective  number  of 
parameters  us  reduced  from  four  to  three. 

Comparison  with  earlier  bandwidth  estimates 

The  ERBs  of  the  filters  obtained  in  the  present  study 
are  compared  with  results  from  three  previous  studies  in 
Figure  4-10  (Zwicker  & Terhardt,  1980;  Moore  & Glasberg,  1983; 
Glasberg  & Moore,  1990)  . The  dotted  lines  are  predictions 
from  the  three  studies.  The  equations  for  the  prediction  of 
critical  bandwidth  are  based  on  predominately  low-frequency 
data.  The  extrapolated  regions  are  shown  as  asterisks.  There 
are  large  difference  in  the  estimated  values  of  the  critical 
bandwidth  obtained  from  these  three  studies,  especially  in  the 
extrapolated  regions.  For  the  present  study,  the  mean  ERB 
values  at  40  dB  SPL  are  shown  as  solid  circles.  The  ERB  value 
at  4 kHz  is  close  to  the  results  of  Moore  and  Glasberg  (1983) 
and  Glasberg  and  Moore  (1990)  . It  is  smaller  than  that  of 
Zwicker  and  Terhardt  (1980)  . This  may  be  due  to  the  different 
methods  used  to  estimate  the  critical  bandwidth.  The  mean  ERB 
values  at  12  and  14  kHz  are  consistent  with  the  prediction  of 
the  modified  formula  given  by  Glasberg  and  Moore  (1990)  . 
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Figure  4-10.  Comparison  of  the  equivalent  rectangular 
bandwidth  (ERB)  of  the  auditory  filter  obtained  in  this  study 
and  the  results  from  three  previous  studies.  The  dotted  lines 
are  given  by  the  formulas  in  the  three  studies.  The  asterisks 
are  the  extrapolated  portions  of  the  formulas.  The  solid 
circles  show  the  mean  ERBs  for  the  three  subjects  in  this 
study. 
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Conclusions 

(1)  Similar  to  the  data  at  low  center  frequencies,  the 
ERB  of  the  auditory  filter  at  12  and  14  kHz  increases  with 
noise  spectrum  level . 

(2)  The  ERBs  of  the  auditory  filter  at  12  and  14  kHz  are 
close  to  the  extrapolated  values  given  by  Glasberg  and  Moore 
(1990)  . 

Experiment  2 

Introduction 

Helmholtz  (1895)  claimed  that  the  ear  is  deaf  to  the 
phase  of  a complex  tone.  He  stated  that  "...  the  quality  of 
the  musical  portion  of  a compound  tone  depends  solely  on  the 
number  and  relative  strength  of  its  partial  simple  tones,  and 
in  no  respect  no  their  difference  of  phase." 

Later  experiments  (Mathes  & Miller,  1947;  Zwicker,  1952; 
Feldtkeller  & Zwicker,  1956;  Goldstein,  1967)  have  shown 
exceptions  to  Helmholtz's  generalization.  The  ear  is 
sensitive  to  the  relative  phase  of  the  components  in  a complex 
tone  if  all  components  of  a complex  fall  within  a single 
auditory  filter. 

The  envelope  detector  model  proposed  by  Goldstein  (1967) 
provides  a good  way  to  understand  the  relationship  between  the 
monaural  phase  discrimination  and  the  auditory  filters.  He 
investigated  two  three -tone  complexes  with  the  same  power 
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spectra,  but  different  relative  phase  among  the  components. 
If  all  three  components  fall  in  the  same  critical  band,  the 
interaction  between  the  components  may  cause  a discernible 
difference  in  envelope  or  time  structure  for  the  two 
complexes.  If  all  the  components  fall  in  different  critical 
band,  the  phase  spectra  of  the  components  cannot  be  perceived 
because  each  tone  is  resolved  by  separate  auditory  filter,  and 
it  is  impossible  to  discriminate  among  the  two  complexes. 

Two  types  of  signals  used  by  Goldstein  (1967)  were 
employed  in  this  study.  One  is  an  amplitude -modulated  (AM) 
signal . It  may  be  expressed  as 

sam=  (l+mcos(umt)  ) cos  (co  ct) 

= {m/2)  cos  ( co  c— to m)  t+cos  (co  ct)  + {m/2)  cos  (to  c+co  m)  t 

Eq.  4-5 

where  <oc  and  com  are  the  angular  frequency  of  the  carrier  tone 
and  the  modulation  tone  respectively,  and  m is  the  modulation 
index.  Another  is  a quasi-frequency-modulated  (QFM)  signal, 
which  may  be  expressed  as 

SQFM=  {m/2)  cos  (co  c-com)  t+cos  (co  ct-7t/2)  + (m/2)  cos  (co  c+com)  t 

Eq.  4-6 

The  AM  and  QFM  signal  have  exactly  the  same  spectral 
components.  The  only  difference  is  a 7r/2  shift  in  the  phase 
of  the  carrier  component.  Figure  4-11  demonstrates  the 
waveforms  and  spectra  of  these  two  signals.  The  envelope 
fluctuation  for  the  AM  signal  (top  panel)  is  much  larger  than 
that  for  the  QFM  signal  (bottom  panel) . In  this  experiment, 
the  subjects  are  required  to  discriminate  the  AM  from  the  QFM 
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signal  for  different  modulation  frequencies.  For  small 
modulation  frequencies,  all  three  components  will  fall  within 
the  same  critical  band  and  the  listener  can  discriminate  the 
envelope  differences  apparent  in  Fig.  4-11.  For  large 
modulation  frequency,  the  three  components  will  fall  into 
different  critical  bands  and  the  relative  phase  of  the 
components  is  inaudible.  The  frequency  boundary  between  the 
two  regions  provides  an  estimate  of  the  critical  band,  or 
frequency  resolving  properties  of  the  ear. 

Method 

Subi ects 

Three  subjects  with  normal  hearing  sensitivity  (250  to 
8000  Hz)  participated  in  this  experiment.  Subject  S3  had 
participated  in  the  previous  experiment.  Subjects  SI  and  S2 
were  new  to  the  experiment  and  received  several  hours  of 
training  in  the  task  before  data  were  collected.  The  average 
age  of  the  subjects  was  20.  Subjects  were  paid  an  hourly  wage 
and  received  a 20  percent  bonus  upon  completion  of  the  study. 
Stimuli 

The  AM  and  QFM  stimuli  were  generated  digitally  using  an 
array  processor  (TDT  QAP1)  controlled  by  an  IBM  PC  compatible 
computer.  The  digital  signal  was  delivered  via  a digital -to- 
analog  converter  (TDT  QDA2)  to  a programmable  attenuator  and 
a digital  switch.  The  attenuated  signal  was  coupled  to  the 
ear  canal  through  the  sound  deliver  system  described  in 
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WAVEFORM 


SPECTRUM 


Figure  4-11.  The  waveform  and  spectrum  of  fully  modulated  (m 
= 1)  AM  and  QFM  signals. 
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Chapter  2.  In  order  to  compensate  for  the  attenuation 
resulting  from  the  sound  delivery  system,  the  input  electrical 
signal  to  the  transducer  was  inversely  filtered  by  the 
frequency  response  of  the  sound  delivery  system,  as  described 
in  the  last  experiment . 

Fully  modulated  (m  = 1)  AM  and  QFM  stimuli  were  employed 
having  carrier  frequencies  of  1,  2,  4,  8,  12,  and  16  kHz.  The 
stimulus  had  a duration  of  500  ms  with  a 20  ms  rise/fall  time. 
The  overall  level  of  the  stimulus  was  chosen  to  be  20  dB 
sensational  level  to  avoid  distortion  products.  This 
precaution  will  be  discussed  later. 

Procedure 

When  collecting  pilot  data,  it  was  found  that  the 
perceptual  basis  for  discriminating  the  AM  from  the  QFM  signal 
varied  depending  on  the  carrier  and  modulation  frequency. 
Thus  a third,  'cue'  interval,  was  added  to  the  forced-choice 
task  to  help  the  subjects  to  distinguish  the  AM  from  the  QFM 
signal.  In  the  modified  forced-choice  procedure,  a trial 
consisted  of  a cue  interval  followed  by  two  choice  intervals. 
The  AM  signal  was  always  presented  in  the  cue  interval.  The 
AM  and  QFM  stimuli  were  presented  in  the  two  intervals  that 
followed  in  a random  order.  The  task  was  to  choose  whether 
the  second  or  third  interval  had  the  AM  signal.  The  between- 
interval  duration  was  300  ms.  A block  of  60  trials  always 
began  with  a small  modulation  frequency.  After  three 
consecutive  correct  responses,  the  modulation  frequency 
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increased  and  for  each  incorrect  response,  the  modulation 
frequency  decreased.  The  step  size,  S,  was  specified  on  a 
logarithmic  scale,  and  expressed  as  101og10  (Afm/fra)  , where  fm  is 
the  current  modulation  frequency,  and  Afra  is  the  change  of 
modulation  frequency.  At  the  beginning  of  a block,  the  step 
size  S was  3.  After  two  reversals,  it  was  decreased  to  1. 
This  three-down  one-up  procedure  estimates  the  79.4%  point  on 
the  psychometric  function  (Levitt,  1971).  Each  critical 
modulation  frequency  point  was  the  average  of  8 blocks. 

As  described  in  Chapter  2,  precautions  were  taken  to 
insure  that  the  sound  source  was  at  the  same  relative  position 
to  the  eardrum  when  the  measurement  was  carried  out . A 
frequency  response  of  the  ear  canal  was  taken  at  the  beginning 
and  end  of  each  session.  For  all  the  data  reported,  the  two 
responses  had  an  rms  difference  of  less  than  1.5  dB  in  the 
frequency  range  from  5 to  16  kHz. 

Results 

The  critical  modulation  frequency  for  the  discrimination 
of  the  AM  and  QFM  signal  is  shown  in  Fig.  4-12.  The  error 
bars  in  the  individual  panels  represent  plus  minus  one 
standard  error  across  8 blocks  of  measurement . The  error  bars 
in  the  average  panel  represent  the  standard  errors  across  the 
means  of  the  three  subjects.  There  are  large  individual 
differences  across  the  three  subjects.  The  average  critical 
modulation  frequency  clearly  increases  with  frequency  from  1 
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to  8 kHz,  but  the  rate  of  increases  diminishes  above  8 kHz. 
In  addition  to  our  estimates,  the  dotted  lines  are  estimates 
of  the  critical  bandwidth  based  on  a variety  of  different 
studies,  including  phase  perception,  loudness  growth  and 
masking  (Zwicker  & Terhardt , 1980;  Moore  & Glasberg,  1983; 
Glasberg  & Moore,  1990) . The  equations  for  the  prediction  of 
critical  bandwidth  are  based  on  predominately  low-frequency 
data.  The  extrapolated  regions  are  shown  as  asterisks.  There 
are  large  difference  in  the  estimated  values  of  the  critical 
bandwidth  obtained  from  these  three  studies,  especially  in  the 
extrapolated  regions.  This  may  due  to  the  different  methods 
they  used  to  estimate  the  critical  band.  Individual 
differences  among  subjects  may  also  contribute  to  the 
discrepancies  in  the  different  estimates  of  critical 
bandwidth.  For  frequencies  below  8 kHz,  the  average  critical 
modulation  frequency  for  our  three  subjects  is  similar  to  the 
critical  bandwidth  curves  estimated  from  the  three  studies. 
Above  8 kHz,  however,  our  average  results  are  much  smaller 
than  the  extrapolated  predictions  of  the  critical  bandwidth. 

Discussion 


Comparison  with  results  from  previous  studies 

Using  a categorical  procedure  somewhat  different  from 
the  one  used  in  our  experiment.  Goldstein  (1967)  obtained  the 
boundaries  for  AM/QFM  signal  discrimination.  His  results 
showed  that  the  boundary  for  AM/QFM  discrimination  was 
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CARRIER  FREQUENCY,  Fc  (kHz) 


Figure  4-12.  The  critical  modulation  frequency  for  the  AM  and 
QFM  signal  discrimination.  The  individual  and  mean  data  are 
shown  as  solid  lines.  The  error  bar  in  the  individual  panel 
represents  plus  minus  one  standard  error  across  eight  blocks 
of  measurement.  The  error  bar  in  the  "mean"  panel  represents 
the  standard  deviation  across  the  means  of  the  three  subjects. 
The  critical  bandwidths  given  by  three  previous  studies  were 
plotted  as  dotted  lines  in  the  "mean"  panel.  The  asterisks 
are  the  extrapolated  portions  of  the  dotted  lines. 
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proportion  to  critical  bandwidth  for  carrier  frequencies  under 
8 kHz.  For  carrier  frequencies  above  8 kHz,  the  boundary  was 
smaller  than  the  extrapolated  prediction  of  the  critical 
bandwidth  estimated  by  Zwicker  and  Terharht  (1980) . This  is 
similar  to  the  results  in  the  present  experiment. 

At  very  small  modulation  index,  the  QFM  signal  is  a very 
good  approximation  of  the  FM  signal.  Feldtkeller  and  Zwicker 
(1956)  studied  the  perception  of  AM  and  FM  signals,  we  can, 
therefore,  their  results  with  the  results  obtained  in  this 
study.  Feldtkeller  and  Zwicker  (1956)  used  a different 
experimental  procedure  to  determine  the  boundaries  of  phase 
perception  as  a function  of  carrier  frequency.  They  estimated 
the  critical  modulation  frequency  for  1 and  4 kHz  carrier 
frequencies  to  be  approximately  90  and  400  Hz  respectively. 
The  critical  modulation  frequency  in  this  study  is 
approximately  125  and  675  Hz  for  1 and  4 kHz  carrier 
frequencies.  Feldtkeller  and  Zwicker' s estimates  are  about 
40%  lower  at  1000  Hz,  and  about  70%  lower  at  4000  Hz  than  our 
estimates.  The  discrepancy  may  be  the  result  of  different 
procedures,  difference  in  listeners,  or  the  influence  of 
distortion  products  as  discussed  in  the  following  section. 

Combination  tones  and  AM/OFM  discrimination 

Later  research  has  implicated  combination  tones  as 
important  contributors  to  the  discrimination  between  AM  and 
QFM  waveforms  (Buunen  & Bilsen,  1974;  Buunen,  Festen,  Bilsen, 
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& van  den  Brink,  1974;  Buunen,  1975) . Buunen  has  suggested 
that  the  discrimination  between  AM  and  QFM  waveforms  is  based 
upon  changes  in  the  internal  spectrum,  which  is  the  result  of 
interactions  between  internally  generated  combination  tones 
and  externally  generated  acoustic  components. 

Smoorenburg  (1972a, b)  has  shown  that  combination  tones 
can  be  perceived  as  long  as  the  frequency  separation  between 
two  tones  is  not  too  large.  Smoorenburg  also  showed  that  the 
audibility  limits  for  perceiving  combination  tones  depended 
upon  the  overall  level  of  the  primaries.  The  phase  and  level 
of  combination  tones  also  depend  on  individual  listeners 
(Humes,  1980) . The  most  widely  studied  combination  tone  is 
the  cubic  difference  tone  (CDT,  2f1-f2)  . For  AM  and  QFM 
stimuli,  one  CDT  is  produced  by  the  interaction  between  the 
middle  (carrier)  and  the  higher  components  and  produces  a 
distortion  product  having  a frequency  exactly  equal  to  the 
frequency  of  the  lower  component  of  the  complex.  The  CDT  and 
the  externally  generated  lower  component  change  what  is  called 
the  ' internal  spectrum' . As  the  phase  of  the  middle  (carrier) 
component  is  different  for  AM  and  QFM  stimuli,  the  CDT  also 
has  a different  phase  for  them.  The  CDT  interacts  differently 
with  the  lower  component  of  the  complex  and  results  in 
different  internal  spectrum  for  AM  and  QFM  stimuli. 

According  to  Smoorenburg  (1972a) , the  CDT  is  barely 
audible  when  the  level  of  the  primaries  are  below  2 0 dB 
sensational  level.  In  order  to  suppress  the  cues  resulting 
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from  the  combination  tones,  we  chose  stimulus  levels  in  this 

experiment  to  be  20  dB  sensation  level.  Thus,  we  are 

confident  that  our  results  are  not  influenced  by  cubic 

difference  tones.  They  could  have  played  a role  in  the 

results  obtained  by  other  investigations. 

Is  AM/OFM  discrimination  a crood  estimate  of  critical 
bandwidth? 

It  is  well  accepted  that  limited  auditory  frequency 
resolution  is  responsible  for  monaural  phase  discrimination. 
Nevertheless,  a simple  concept  of  frequency  resolution  as 
measured  by  the  critical  band  is  inadequate  to  account  for  all 
of  the  phase  discrimination  results.  For  different  signal 
levels,  Goldstein  (1967)  has  shown  that  the  critical 
modulation  frequency  for  AM/QFM  discrimination  increased  as  a 
function  of  the  signal  level.  This  is  consistent  with 
evidence  that  the  critical  band  widens  as  a function  of  level 
(Pick,  1980;  Lutfi  & Patterson,  1984;  Patterson  & Moore, 
1986)  . Alternatively,  the  influence  of  combination  tones  also 
increase  with  level,  and  they  may,  therefore  change  the 
critical  modulation  frequency. 

For  frequencies  above  8 kHz,  another  factor  may 
contribute  to  the  discrepancy  between  the  measured  critical 
modulation  frequency  and  the  predicted  critical  bandwidth. 
For  these  high  modulation  frequencies,  even  if  all  components 
of  the  complex  fall  within  the  same  critical  band,  the 
auditory  process  may  be  too  sluggish  to  follow  the  resulting 
envelope  patterns.  Viemeister  (1979),  for  example,  measured 
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the  temporal  modulation  transfer  function  (TMTF)  of  a 
sinusoidally  amplitude-modulated  wide  band  noise.  His  results 
showed  that  observers  failed  to  detect  temporal  fluctuations 
for  modulation  frequencies  above  1000  Hz.  Frisina,  Smith  and 
Chamberlain  (1990)  investigated  the  encoding  of  amplitude 
modulation  in  a single  neuron  of  the  auditory  nerve  and 
ventral  cochlear  nucleus  (VCN) , also  showed  that  the  highest 
modulation  frequency  that  could  be  encoded  was  approximately 
500  Hz. 

Conclusion 

The  present  experiment  shows  that  for  low  stimulus  level, 
the  critical  modulation  frequency  for  monaural  phase 
discrimination  is  roughly  in  agreement  with  the  critical 
bandwidth  estimated  by  the  previous  studies  for  frequencies 
below  8 kHz.  For  frequencies  above  8 kHz,  the  critical 
modulation  frequency  is  significantly  smaller  than  the 
predicted  critical  bandwidth.  Although  the  limited  frequency 
resolution  is  responsible  for  the  monaural  phase 
discrimination,  there  are  not  enough  evidence  to  show  that  the 
critical  band  obtained  from  masking  experiments  is  the  sole 
mechanism  for  the  phase  discrimination. 


Chapter  5 
SUMMARY 

The  measurement  of  high-frequency  hearing  is  very 
important  in  the  monitoring  of  various  hearing  losses.  The 
present  study  dealt  with  two  primary  issues  involved  in  high- 
frequency  hearing:  absolute  sensitivity  and  frequency 
resolution . 

The  primary  problem  with  the  measurement  of  high- 
frequency  sensitivity  arises  from  the  fact  that  standing  waves 
are  generated  in  the  ear  canal  for  frequencies  above  3500  Hz. 
Because  of  the  slope  of  the  eardrum,  the  resonances  caused  by 
the  standing  waves  result  in  different  distributions  of  sound 
pressure  over  the  surface  of  the  tympanic  membrane  at 
different  signal  frequencies.  The  purpose  of  experiments  in 
Chapter  2 was  to  determine  the  relationship  between  the 
distribution  of  sound  pressure  and  the  audiometric  absolute 
threshold . 

In  order  to  investigate  the  effect  of  standing  waves  on 
the  estimates  of  high-frequency  thresholds,  an  audiometric 
system  was  designed  to  transmit  sound  from  a transducer  to  the 
ear  canal  by  a long  plastic  tube  with  a diameter  equal  to  that 
of  the  average  ear  canal.  In  one  experimental  condition,  the 
long  plastic  tube  was  used  as  the  sound  source  to  drive  the 
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ear  canal.  Because  of  the  loss  in  the  long  plastic  tube,  this 
system  produced  little  reflection  at  the  entrance  to  the  ear 
canal . Thus  there  was  essentially  no  change  in  sound  pressure 
as  a function  of  frequency  at  the  position  of  the  eardrum  with 
a fixed  input  level.  In  another  experimental  condition,  a 
hard-wall  material  was  inserted  in  the  plastic  tube  at  a 
position  outside  the  entrance  of  the  ear  canal.  This  created 
large  changes  in  sound  pressure  at  the  eardrum.  The 
difference  in  threshold  obtained  in  these  two  different 
conditions  revealed  the  effect  of  change  in  sound  pressure  at 
the  eardrum  resulting  from  standing  waves  in  the  ear  canal. 
The  results  showed  that,  in  the  frequency  range  from  5 to  16 
kHz,  the  change  in  absolute  threshold  as  a function  of 
frequency  caused  by  the  standing  waves  was  about  ± 5.0  dB . 

A physical  model  was  proposed  to  predict  the  sound 
pressure  at  the  eardrum  from  the  measured  sound  pressure  at 
the  entrance  of  the  ear  canal.  It  was  shown  that  for  a fixed 
input  level,  a good  estimate  of  the  change  in  pressure  as  a 
function  of  frequency  at  the  eardrum  (termed  the  MMSPR)  was 
one-half  the  MMSPR  measured  at  the  entrance  of  the  ear  canal 
on  a logarithmic  scale.  The  results  showed  that  the  predicted 
sound  pressure  at  the  eardrum  correlated  well  with  the 
changing  pattern  of  the  absolute  threshold  at  frequencies 
below  13  kHz.  It  was  concluded  that  a long  lossy  tube 
provided  a simple  means  of  minimizing  the  changes  in  pressure 
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at  the  eardrum  created  by  strong  resonances  and  anti- 
resonances  in  the  ear  canal . 

Because  standing  waves  result  in  different  sound-pressure 
levels  at  the  eardrum  for  different  positioning  of  the 
transducer,  this  produces  variations  between  the  test-retest 
results.  In  Chapter  3,  the  fitting  variance  of  intra-subject 
threshold  caused  by  transducer  placement  was  estimated  for  two 
types  of  transducers,  the  hard-wall  source  and  headphones. 
For  the  hard-wall  source,  the  variation  in  behavior  threshold 
was  0.7  dB  for  a 1-kHz  signal  and  2.9  dB  for  signals  above  7 
kHz  if  the  sound  source  moved  by  ± 0.5  cm.  The  smaller 
variation  in  low-frequency  threshold  was  due  to  its  larger 
wavelength.  It  was  found  out  that  the  fitting  variance  for 
high-frequency  signals  was  much  larger  than  that  of  the  low- 
frequency  signals. 

For  headphones,  it  was  found  that  the  sound  pressure 
change  in  the  ear  canal  caused  by  headphone  placement  was 
smaller  than  that  caused  by  hard-wall  source  movement.  This 
was  probably  due  to  the  smaller  acoustic  impedance  for  the 
headphones.  Thus  for  high-frequency  signals,  the  resulting 
fitting  variance  measured  with  the  headphones  was  expected  to 
be  smaller  than  + 2.9  dB  as  measured  with  high-wall  source. 

The  effect  of  the  slope  of  the  psychometric  function  on 
the  threshold  variance  was  also  estimated.  It  was  found  that 
the  slope  of  the  psychometric  function  for  a 15 -kHz  signal  was 
25%  larger  than  that  of  a 1-kHz  signal.  This  may  result  in 
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25%  smaller  threshold  variance  for  high-frequency  than  low- 
frequency  signals.  The  larger  fitting  variance  and  smaller 
threshold  variance  for  high-frequency  signals  may  be 
responsible  for  the  diverse  reports  on  the  reliability  of 
high-frequency  threshold  estimates. 

In  Chapter  4,  we  estimated  the  frequency  resolution  of 
the  ear  for  frequencies  above  10  kHz  using  two  previously 
established  methods.  The  first  method  employed  notched-noise 
maskers.  The  results  showed  that  at  12  and  14  kHz,  the 
bandwidths  of  the  estimated  auditory  filters  increased  with 
noise  spectrum  level.  This  was  consistent  with  the  results 
obtained  for  lower  center  frequencies.  The  equivalent 
rectangular  bandwidths  of  the  estimated  auditory  filters  at  12 
and  14  kHz  were  close  to  the  extrapolated  values  given  by 
Glasberg  and  Moore  (1990) . The  second  method  was  based  on  the 
monaural  phase  discrimination  of  amplitude -modulated  (AM)  and 
quasi- frequency-modulated  (QFM)  signals.  For  low  stimulus 
level,  the  measured  critical  modulation  frequency  at  center 
frequencies  below  8 kHz  was  roughly  in  agreement  with  the 
critical  bandwidth  obtained  with  the  masking  method.  At 
center  frequencies  above  8 kHz,  however,  the  measured  critical 
modulation  frequency  was  smaller  than  the  critical  bandwidth 
obtained  with  the  masking  method.  According  to  the  envelope 
detecting  model  proposed  by  Goldstein  (1967) , the 
discrimination  of  AM  and  QFM  signals  was  based  on  the  envelope 
fluctuation  pattern  when  all  the  components  fall  in  the  same 
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critical  band.  For  very  high  center  frequencies,  the  auditory 
process  may  be  too  sluggish  to  follow  the  envelope 
fluctuations  even  when  all  components  fall  into  the  same 
critical  band.  This  may  contribute  to  the  discrepancy  between 
the  results  obtained  with  the  phase  discrimination  and  the 
masking  method  at  very  high  center  frequencies . 


REFERENCES 


Barr-Hamilton,  R.M. , Bryan,  M.E.,  & Tempest,  W.  (1968).  The 

application  of  signal  detection  theory  to  the  quiet 
threshold  of  hearing.  Sound.  2,  92-98. 

Butler,  R.A.,  & Helwig,  C..C.  (1983).  The  spatial  attributes 

of  stimulus  frequency  in  the  median  sagittal  plane  and 
their  role  in  sound  localization.  American  Journal  of 
Otolaryngology . 4,  165-173. 

Buunen,  T.J.F.  (1975) . Two  hypotheses  on  monaural  phase 
effects.  Acustica . 34 , 98-105. 

Buunen,  T.J.F. , & Bilsen,  F . A . (1974).  Subjective  phase 

effects  & combination  tones,  in  Facts  and  Models  in 
Hearing,  edited  by  E.  Zwicker  and  E.  Terhardt  (Springer- 
Verlag,  New  York) . 

Buunen,  T.J.F.,  Festen,  J.M.,  Bilsen,  F . A . & van  den  Brink,  G. 
(1974).  Phase  effects  in  a three -component  signal. 
Journal  of  the  Acoustical  Society  of  America,  55 , 297- 

303  . 

Carlyon,  R.P.  (1989).  Changes  in  the  masked  thresholds  of 
brief  tones  produced  by  prior  bursts  of  noise.  Hearing 
Research . 41,  223-236. 

Collins,  M.,  Cullen,  J. , & Berlin  C.  (1981).  Auditory  signal 
processing  in  a hearing  impaired  subject  with  residual 
ultra-audiometric  hearing.  Audiology,  20 , 347-361. 

Dai,  H.  (1994).  Signal-frequency  uncertainty  in  spectral-shape 
discrimination:  Psychometric  functions.  Journal  of  the 

Acoustical  Society  of  America  (in  review) . 

Dallos,  P.  (1973).  The  auditory  periphery.  Biophysics  and 
physiology.  New  York:  Academic  Press. 

De  Seta,  E.,  Bertoli,  G.A. , & Filipo,  R.  (1985).  High- 

frequency  audiometry  above  8 kHz.  Audiology,  24 , 254-259. 

Dreisbach,  L.E.,  & Frank,  T.  (1991).  Low-  and  high-frequency 
threshold  reliability.  Journal  of  the  Acoustical  Society 
of  America.  89 . 1975. 


132 


133 


Dreschler,  VI. A.,  van  der  Hulst,  R.J.A.M.,  Tange,  R.A.,  & 

Urbanus,  N.A.M.  (1985) . The  role  of  high-frequency 
audiometry  in  early  detection  of  ototoxicity.  Audiology. 
24 . 387-395. 

Duifhuis,  H.  (1970) . Audibility  of  high  harmonics  in  a 
periodic  pulse.  Journal  of  the  Acoustical  Society  of 
America . 48 . 888-893. 

Duifhuis,  H.  (1971) . Audibility  of  high  harmonics  in  a 
periodic  pulse.  II,  Time  effect.  Journal  of  the 
Acoustical  Society  of  America,  49.,  1155-1162. 

Egan  J.P.,  Lindner,  VI. A.,  & McFadden,  D.  (1966).  Masking-level 
differences  and  the  form  of  the  psychometric  function. 
Journal  of  the  Acoustical  Society  of  America.  37., 
1181 (A) . 

Elliot,  P.B.  (1964).  Tables  of  d'  . in  Signal  detection  and 
recognition  by  human  observers,  J.A.  Swets  (Ed.),  651-684 
(New  York,  Wiley) . 

Fausti,  S.A.,  Erickson,  D.A.,  Frey,  R.H.,  Rappaport,  B.Z.,  & 
Schechter,  M . A . (1981) . The  effects  of  noise  upon  human 

hearing  sensitivity  from  8000-20000  Hz.  Journal  of  the 
Acoustical  Society  of  America,  6 9 , 1343-1349. 

Fausti,  S.A.,  Frey,  R.H.,  Erickson,  D.A.,  Rappaport,  B.Z., 
Cleary,  E.J.,  & Brummett,  R.E.  (1979).  A system  for 
evaluating  auditory  function  from  8000-20000  Hz.  Journal 
of  the  Acoustical  Society  of  America.  66 , 1713-1718. 

Fausti,  S.A.,  Schechter,  M.A.,  Rappaport,  B.Z.,  Frey,  R.H.,  & 
Mass,  R.E.  (1984) . Early  detection  of  cis-platin 

ototoxicity.  Cancer , 53 , 224-231. 

Feldtkeller,  R.  , & Zwicker,  E.  (1956).  Das  ohr  als 

nachrichtenempf anger . (S.  Hirzel  Verlag,  Stuttgart) , 

chap . 4 . 

Fletcher,  H.  (1940).  Auditory  patterns.  Review  of  Modern 
Physics , 12.,  47-65. 

Fletcher,  J.L.  (1965).  Reliability  of  high-frequency 

thresholds.  Journal  of  Auditory  Research.  5.,  133-137. 

Fletcher,  J.L.,  Cairns,  A.B.,  Collins,  F.G.,  & Endicott,  J. 

(1967) . High  frequency  hearing  following  meningitis. 
Journal  of  Auditory  Research,  1_,  223-227 . 


134 


Foster,  S.  (1986).  Impulse  response  measurement  using  Golay 
codes.  IEEE,  Acoustics.  Speech  and  Signal  Processing, 
929-932 . 

Frisina  R.D.,  Smith,  R.L.,  & Chamberlain,  S.C.  (1990). 

Encoding  of  amplitude  modulation  in  gerbil  cochlear 
nucleus:  II  possible  neural  mechanisms.  Hearing  Research, 
44 , 123-142. 

Gardner,  M.B.,  & Gardner,  R.S.  (1973).  Problem  of  localization 
in  the  medial  plane:  effect  of  pinnae  cavity  occlusion. 
Journal  of  the  Acoustical  Society  of  America.  53 . 400- 

408  . 

Glasberg,  B.R.,  & Moore,  B.C.J.  (1986).  Auditory  filter  shapes 
in  subjects  with  unilateral  and  bilateral  cochlear 
impairments.  Journal  of  the  Acoustical  Society  of 
America , 79 . 1020-1033. 

Glasberg,  B.R.,  & Moore,  B.C.J.  (1990).  Derivation  of  auditory 
filter  shapes  from  notched-noise  data.  Hearing  Research, 
47 , 103-138. 

Glasberg,  B.R.,  Moore,  B.C.J.,  Patterson,  R.D.,  & Nimmo-Smith, 
I.  (1984).  Dynamic  range  and  asymmetry  of  the  auditory 
filter.  Journal  of  the  Acoustical  Society  of  America,  76, 
419-427 . 

Golay,  M.J.E.  (1961) . Complementary  series.  IRE  Transactions 
on  Information  Theory,  1_,  82-87. 

Goldstein,  J.L.  (1967) . Auditory  spectral  filtering  and 
monaural  phase  perception.  Journal  of  the  Acoustical 
Society  of  America.  41 . 458-479. 

Green,  D.M.  (1990) . Stimulus  selection  in  the  adaptive 
psychophysical  procedures.  Journal  of  the  Acoustical 
Society  of  America,  87 , 2662-2667. 

Green,  D.M.,  & Dai,  H.  (1991).  Probability  of  being  correct 

with  1 of  M orthogonal  signals.  Perception  and 
Psychophysics . 49 , 100-101. 

Green,  D.M,  Kidd,  G.  , Jr.,  & Stevens,  K.H.  (1987).  High- 
frequency  audiometric  assessment  of  a young  adult 
population.  Journal  of  the  Acoustical  Society  of  America, 
81 . 485-494. 


Green  D.M. , & Swets,  J.A.  (1966)  . Signal  detection  theory  and 
psychophysics . New  York:  John  Wiley  & Sons. 


135 


Hamilton,  P.M.  (1957) . Noise  masked  thresholds  as  a function 
of  tonal  duration  and  masking  noise  band  width.  Journal 
of  the  Acoustical  Society  of  America,  2 9 . 506-511. 

Harris,  J.D.,  & Myers,  C.K.  (1971).  Tentative  audiometric 

threshold-level  standards  from  8 through  18  kHz.  Journal 
of  the  Acoustical  Society  of  America,  49 , 600-601. 

Helmholtz,  H.L.F.  (1895) . Sensation  of  tone,  pp  126-127 

Henry,  K.R.,  Fast,  G.A.,  Nguyen,  H.H.,  Paolinelli,  M.C.,  & 

Ayars,  N.M.  (1985) . Extra-high- frequency  auditory 

thresholds:  Fine  structure,  reliability,  temporal 

integration  and  relation  to  the  ear  canal  resonance. 
Audiology.  24 , 92-103. 

Houtgast,  T.  (1974) . Lateral  suppression  in  hearing,  doctoral 
dissertation.  (University  of  Amsterdam,  Amsterdam,  The 
Netherlands) . 

Hudde,  H.  (1983) . Measurement  of  the  eardrum  impedance  of 
human  ears . Journal  of  the  Acoustical  Society  of  America, 
73  . 242-247. 

Humes,  L.E.  (1980) . On  the  nature  of  tow- tone  aural 
nonlinearity.  Journal  of  the  Acoustical  Society  of 
America . 67 , 2073-2083. 

Jacobson,  E.J.,  Downs,  M.P.,  & Fletcher,  J.L.  (1969).  Clinical 
findings  in  high-frequency  thresholds  during  known 
ototoxic  drug  usage.  Journal  of  Auditory  Research,  9, 
379-385 . 

Johansen,  P.A.  (1975).  Measurement  of  the  human  ear  canal. 
Acustica , 33 , 349-351. 

Khanna,  S.M.,  & Stinson,  M.R.  (1985).  Specification  of  the 

acoustical  input  to  the  ear  at  high  frequencies.  Journal 
of  the  Acoustical  Society  of  America.  77,  577-589. 

Killion,  M.C.  (1978).  Revised  estimate  of  minimum  audible 
pressure.  Where  is  the  'missing  6 dB  ?' . Journal  of  the 
Acoustical  Society  of  America.  63 , 1501-1508. 

Leek,  M.R.,  Hanna,  T.E.,  & Marshall,  L.  (1992).  Estimation  of 
psychometric  functions  from  adaptive  tracking  procedures. 
Perception  and  Psychophysics,  51 , 247-256. 

Levitt,  H.  (1971).  Transformed  up-down  methods  in 
psychoacoustics.  Journal  of  the  Acoustical  Society  of 
America , 49 . 467-477. 


136 


Lutfi,  R.A.,  & Patterson,  R.D.  (1984).  On  the  growth  of 

masking  asymmetry  with  stimulus  intensity.  Journal  of  the 
Acoustical  Society  of  America,  76 , 739-745. 

Mathes,  R.C.,  & Miller,  R.L.  (1947).  Phase  effects  in  monaural 
phase  perception.  Journal  of  the  Acoustical  Society  of 
America , 19 . 780-797. 

Middlebrooks , J.  C.,  & Green,  D.  M.  (1990).  Directional 
dependence  of  interaural  envelope  delays.  Journal  of  the 
Acoustical  Society  of  America.  87.,  2149-2162. 

Moore,  B.C.J.,  & Glasberg,  B.R.  (1983)  Suggested  formulae  for 
calculating  auditory-filter  bandwidths  and  excitation 
patterns.  Journal  of  the  Acoustical  Society  of  America, 
74,  750-753. 

Neumann  J.  , Uppenkamp,  S.,  & Kollmeier,  B.  (1993). 

Chirpsignals  with  low  peakfactors  and  freely  selectable 
power  spectra.  Contributions  to  Psychological  Acoustics, 
Bibliotheks-  und  Inf ormat ionssystem  der  Universitat 
Oldenburg,  89-95. 

Northern,  J.L.,  Downs,  M.P.,  Rudmose,  W.  , Glorig,  A.,  & 

Fletcher,  J.L.  (1972) . Recommended  high-frequency 
audiometric  threshold  levels  (8000-18000  Hz) . Journal  of 
the  Acoustical  Society  of  America.  52.,  585-595. 

Ohm,  G.S.  (1843).  Uber  die  definition  des  tones,  nebst  daran 
gekniipfter  theorie  der  sirene  und  ahnlicher  tonbildender 
vorrichtungen . Annals  of  Physical  Chemistry,  59.,  513-565. 

Osterhammel,  D.,  Osterhammel,  P.,  & Terkildsen,  K.  (1977).  A 
quasi-f ree-f ield  transducer  system  for  high-frequency 
audiometry.  Scandinavian  Audiology  6.,  91-95. 

Patterson,  R.D.  (1976) . Auditory  filter  shapes  derived  with 
noise  stimuli.  Journal  of  the  Acoustical  Society  of 
America , 59 , 640-654. 

Patterson,  R.D.,  & Henning,  G.B.  (1977).  Stimulus  variability 
and  auditory  filter  shape.  Journal  of  the  Acoustical 
Society  of  America,  62.,  649-664. 

Patterson,  R.D.,  & Moore,  B.C.J.  (1986).  Auditory  filters  and 
excitation  patterns  as  representations  of  frequency 
resolution,  in  Frequency  selectivity  in  hearing.  London: 
Academic  Press. 

Patterson,  R.D.,  & Nimmo-Smith,  I.  (1980).  Offer  frequency 

listening  and  auditory-filter  asymmetry.  Journal  of  the 
Acoustical  Society  of  America,  67.,  229-245. 


137 


Patterson,  R.D.,  Nimmo-Smith,  I.,  Weber,  D.L.,  & Milroy,  R. 

(1982)  . The  deterioration  of  hearing  with  age:  Frequency 
selectivity,  the  critical  ratio,  the  audiogram,  and 
speech  threshold.  Journal  of  the  Acoustical  Society  of 
America , 72 , 1788-1803. 

Pick,  G.F.  (1980).  Level  dependence  of  psychophysical 
frequency  resolution  and  auditory  filter  shape.  Journal 
of  the  Acoustical  Society  of  America,  68 , 1085-1095. 

Plomp  R.  (1964) . The  ear  as  a frequency  analyzer.  Journal  of 
the  Acoustical  Society  of  America.  3 6 , 1628-1636. 

Pumplin,  J.  (1985)  Low-noise  noise.  Journal  of  the  Acoustical 
Society  of  America,  78 , 100-104. 

Robbins,  H.,  & Monro,  S.  (1951).  A stochastic  approximation 

method.  Annals  of  Mathematical  Statistics,  22 , 400-407. 

Roffler,  S.K.,  & Butler,  R . A . (1968).  Factors  that  influence 
the  localization  of  sound  in  the  vertical  plane.  Journal 
of  the  Acoustical  Society  of  America,  43.,  1255-1259. 

Rosen,  S.,  Plester,  D,  El-Moffy,  A.,  & Rosen,  H.  (1964).  High- 
frequency  audiometry  in  presbycusis.  Archives  of 
Otolaryngology , 79 , 18-32. 

Scholl,  H.  (1962) . Das  dynamische  verhalten  des  Gehors  bei  der 
Unterteilung  des  schallspektrums  Frequenzgruppen . 
Acustica . 12.,  101-107. 

Schroeder,  M.  R.  (1970) . Synthesis  of  low-peak-factor  signals 
and  binaural  sequences  with  low  autocorrelation.  IEEE 
Transaction  on  Information  Theory,  IT- 16 , 85-89. 

Shailer,  M.J.,  Moore,  B.C.J.,  Glasberg,  B.R.,  Watson,  N.  , 
Harris,  S.  (1990).  Auditory  filter  shapes  at  8 and  10 
kHz.  Journal  of  the  Acoustical  Society  of  America,  88 , 
141-148 . 

Shaw,  E.A.G.  (1966) . Ear  canal  pressure  generated  by 
circumaural  and  supraaural  earphones.  Journal  of  the 
Acoustical  Society  of  America,  39 , 471-479. 

Smoorenburg,  G.F.  (1972a).  Audibility  region  of  combination 
tones.  Journal  of  the  Acoustical  Society  of  America,  52, 
603-614 . 

Smoorenburg,  G.F.  (1972b) . Combination  tones  and  their  origin. 
Journal  of  the  Acoustical  Society  of  America,  52.,  615- 


632  . 


138 


Stelmachowicz , P.G.,  Beauchaine,  K . A . , Kalberer,  A.,  Kelly, 
W.J.,  & Jesteadt,  W.  (1989).  High-frequency  audiometry: 
Test  reliability  and  procedure  considerations.  Journal  of 
the  Acoustical  Society  of  America.  85 , 879-887. 

Stelmachowicz,  P.G.,  Beauchaine,  K . A . , Kalberer,  A.  Langer, 
T.,  & Jesteadt,  W.  (1988).  The  reliability  of  auditory 

thresholds  in  the  8-  to  20 -kHz  range  using  a prototype 
audiometer.  Journal  of  the  Acoustical  Society  of  America, 
83 . 1528-1535. 

Stevens,  K.N.,  Berkovitz,  R.  , Kidd,  G.,  Jr.,  & Green,  D.M. 
(1987) . Calibration  of  ear  canals  for  audiometry  at  high 
frequencies.  Journal  of  the  Acoustical  Society  of 
America , 81 . 470-484. 

Stinson,  M.R.  (1985) . The  spatial  distribution  of  sound 
pressure  within  scaled  replicas  of  the  human  ear  canal. 
Journal  of  the  Acoustical  Society  of  America,  78.,  1596- 
1602  . 

Stinson,  M.R.  (1990) . Revision  of  estimates  of  acoustic  energy 
reflectance  at  the  human  eardrum.  Journal  of  the 
Acoustical  Society  of  America,  88 , 1773-1778. 

Swets,  J.A.,  Green,  D.M.,  & Tanner,  W.P.  (1962).  On  the  width 
of  critical  bands.  Journal  of  the  Acoustical  Society  of 
America , 34 . 108-113. 

Tange,  R.A.,  Conijn,  E.A.J.G.,  & van  Zeijl,  L.G.P.M.  (1982a). 
The  cortitoxic  effects  of  cis-platinum  in  the  guinea  pig. 
Archives  of  Otorhinolaryngology.  237 . 17-26. 

Tange,  R.A.,  Conijn,  E.A.J.G.,  van  Zeijl,  L.G.P.M.,  & Huizing, 
E.H.  (1982b).  Pattern  of  gentamicin- induced  cochlear 
degeneration  in  the  guinea  pig-  Archives  of 

Otorhinolaryngology . 236 , 173-184. 

Taylor,  M.M.  (1971) . On  the  efficiency  of  psychophysical 
measurement . Journal  of  the  Acoustical  Society  of 
America , 49 . 505-507. 

van  Zeijl,  L.G.P.M.,  Conijn,  E.A.J.G.,  Rodenburg,  M. , Tange, 
R.A.  & Brocaar,  M.P.  (1984) . Analysis  of  hearing  loss  due 
to  the  cis-diamminedichloroplatinum- II . Archives  of 
Otorhinolaryngology , 239,  255-262. 

Viemeister,  N.F.  (1979).  Temporal  modulation  transfer 

functions  based  upon  modulation  thresholds.  Journal  of 
the  Acoustical  Society  of  America,  66.,  1364-1380. 


139 


Viemeister,  N.F.  (1982) . Intensity  discrimination  at  high 
frequencies  in  the  presence  of  noise.  Journal  of  the 
Acoustical  Society  of  America,  72 , S89. 

Viemeister,  N.F.,  & Bacon,  S.P.  (1983).  Increment  detection, 
intensity  discrimination,  and  tone  decay  at  high 
frequencies.  Journal  of  the  Acoustical  Society  of 
America , 73 , S93 . 

von  Bekesy,  G.  (1960)  . Experiments  in  hearing.  New  York: 
McGraw  Hill. 

von  Helmholtz,  H.L.F.  (1954) . On  the  sensations  of  tone  as  a 
physiological  basis  for  the  theory  of  music.  New  York: 
Dover.  (Original  work  published,  1863) 

Voss,  S.E.,  & Allen,  J.B.  (1994).  Measurement  of  acoustic 

impedance  and  reflectance  in  the  human  ear  canal.  Journal 
of  the  Acoustical  Society  of  America,  95 , 372-384. 

Watson,  C.S.,  Franks,  J.A.,  & Hood,  D.C.  (1967).  Detection  of 
tones  in  the  absence  of  external  masking  noise.  Journal 
of  the  Acoustical  Society  of  America,  42 , 1194 (A) . 

Watson,  C.S.,  Franks,  J.R.,  & Hood,  D.C.  (1972).  Detection  of 
tones  in  the  absence  of  external  masking  noise.  I.  Effect 
of  signal  intensity  and  signal  frequency.  Journal  of  the 
Acoustical  Society  of  America,  52.,  633-643. 

Wegel,  R.L.  (1932) . Physical  data  and  physiology  of  excitation 
of  the  auditory  nerve.  Annals  of  Otology,  Rhinology,  & 
Laryngology,  41 , 740-779. 

Wegel,  R.L.  & Lane,  C.E.  (1924) . The  auditory  masking  of  one 
pure  tone  by  another  and  its  possible  relation  to  the 
dynamics  of  the  inner  ear.  Physics  Review.  23 . 266-285. 

Wetherill,  G.B.  (1966) . Seouential  methods  in  statistics. 
London:  Methuen. 

Zhou,  B.,  Green,  D.M.,  & Middlebrooks , J.C.  (1992). 

Characterization  of  external  ear  impulse  responses  using 
Golay  codes.  Journal  of  the  Acoustical  Society  of 
America . 92.,  1169-1171. 

Zislis , T.,  Sc  Fletcher,  J.L.  (1966).  Relation  of  high- 

frequency  thresholds  to  age  and  sex.  Journal  of  Auditory 
Research , 6,  189-198. 

Zwicker,  E.  (1952)  . Die  grenzen  der  horbarkeit  der 
amplitudenmodulation  und  der  f requenzmodulation  eine 
tones.  Acustica . 2,  Beih.  3,  125-133. 


140 


Zwicker,  E.  (1956)  . Die  elementaren  grundlagen  zur  bestimmung 
der  inf ormationskapazitat  des  gehors . Acustica , 6.,  365- 
381 

Zwicker,  E.,  & Terhardt , E.  (1980).  Analytical  expressions  for 
critical  band  rate  and  critical  bandwidth  as  a function 
of  frequency.  Journal  of  the  Acoustical  Society  of 
America , 68 . 1523-1525. 


BIOGRAPHICAL  SKETCH 


Bin  Zhou,  born  in  Hubei,  P.  R.  China,  on  August  22,  1967. 
He  graduated  from  Dongfanghong  Senior  High  School  in  1983. 
Bin  received  a B.  S.  in  acoustics  in  1987  from  Nanjing 
University  of  P.  R.  China.  After  graduating  from  college,  he 
worked  as  an  assistant  engineer  at  Kesheng  Technology  in  Wuhan 
for  three  years  before  coming  to  the  United  States . Bin 
enrolled  in  the  doctoral  program  in  psychology  at  the 
University  of  Florida  in  August  of  1990,  where  he  studied 
psychoacoustics  under  the  supervision  of  Professor  David  M. 
Green.  Bin  plans  a career  in  academic  research  of  speech  and 
hearing . 


141 


I certify  that  I have  read  this  study -and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

7ft  iMues 

David  M.  Green,  Chairman 
Graduate  Research  Professor 
of  Psychology 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


6>Wr§  i/fvvU L ^ 

Robert  D.  Sorkin 
Professor  of  Psychology 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Professor  of  Psychology 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Fohn  C.  Middlebrooks 
'Associate  Professor 
of  Neuroscience 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Donald  G.  Childers 
Professor  of  Electrical 
Engineering 


This  dissertation  was  submitted  to  the  Graduate  Faculty 
of  the  Department  of  Psychology  in  the  College  of  Liberal  Arts 
and  Sciences  and  to  the  Graduate  School  and  was  accepted  as 
partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy. 


August,  1994 


Dean,  Graduate  School 


